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NON-RAYLEIGH  SEA  CLUTTER:  PROPERTIES  AND  DETECTION  OF  TARGETS 


I.  INTRODUCTION 

In  the  last  several  years,  many  investigators  have  studied  the  scattering  mechanisms 
that  produce  radar  sea  clutter.  With  results  for  scattering  from  slightly  rough  surfaces 
and  composite  surfaces  obtained  by  Rice  [1],  Wright  [2,3],  Valenzuela  [4-6],  Guinard 
and  Daley  [7],  and  others,  the  properties  of  average  radar  backscatter  can  be  modeled  • 
fairly  well. 

The  original  work  on  the  probability  density  p(x)  of  sea  clutter  is  that  of  Goldstein 
[8].  Goldstein  states  that  if  many  scatterers  are  uniformly  distributed  in  an  illuminated 
patch  (the  area  defined  by  the  pulsewidth  and  the  radar  beamwidth,  cr/2  by  R8),  the 
relative  phases  of  the  individual  echoes  will  be  random.  He  adds  that  the  central  limit 
theorem  yields  the  Rayleigh  density  for  envelope-detected  sea  clutter. 

He  notes,  however,  that  if  the  pulsewidth  is  small,  the  assumption  of  uniformly  dis- 
tributed scatterers  does  not  hold.  As  an  example,  he  shows  a photograph  of  an  A-scope 
(Fig.  1)  and  notes  the  “spiky”  appearance.  Furthermore,  he  infers  that  the  radar  is  in 
fact  resolving  the  individual  waves. 

This  work  was  published  in  1951.  During  the  next  15  years,  very  little  research  was 
done  on  the  density  of  sea  clutter,  and  that  which  was  done  was  either  classified  or  ap- 
peared in  reports  of  limited  circulation.  In  1969,  Nathanson  [9]  reported  some  results 
obtained  by  the  Naval  Research  Laboratory  (NRL)  and  the  Applied  Physics  Laboratory  of 
the  Johns  Hopkins  University  (APL).  Specifically,  he  gave  standard-deviation-to-mean 
ratios  for  various  pulsewidths  and  showed  a deviation  from  the  Rayleigh  density  for  small 
pulse  widths.  Furthermore,  he  stated  that  for  short  pulses,  the  density  function  for  hori- 
zontal polarization  had  a longer  tail  than  the  density  function  for  vertical  polarization. 

In  1970,  Trunk  and  George  [10]  considered  the  log-normal  and  contaminated-normal 
descriptions  of  sea  clutter  and  calculated  detection  probabilities  for  targets  in  these  den- 
sities. With  that,  the  detailed  description  of  the  clutter  density  will  begin.  In  Sec.  II, 
measurements  of  the  average  backscatter  Oq  are  given,  and  the  results  of  slightly  rough 
scattering  and  the  composite  surface-scattering  model  are  introduced  so  that  they  may  be 
used  later  in  this  report.  In  Sec.  Ill,  a spatially  varying  conditional  density  p(x|0o)  is 
introduced.  This  conditional  density  is  a natural  consequence  of  the  composite  surface- 
scattering model  and  is  used  to  explain  the  non-Rayleigh  nature  of  sea  clutter.  The  varia- 
tion of  the  clutter  densities,  p(x)  and  p(x|cro),  with  various  radar  parameters  such  as  fre- 
quency, pulsewidth,  and  polarization  are  found  using  analysis  of  variance  techniques.  In 
Sec.  IV,  a way  of  constructing  a realistic  computer  model  of  the  sea  surface  is  given.  The 
surface  is  used  to  predict  non-Rayleigh  clutter  densities  for  various  conditions  and  to  indi- 
cate some  finer  points  associated  with  the  detection  of  small  targets  on  the  surface  of  the 
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ocean.  In  Sec.  V,  various  detectors  such  as  integrators,  rank  detectors,  and  the  Neyman- 
Pear  a detector  are  compared.  The  major  results  of  the  report  are  summarized  in  Sec.  VI. 


I 
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II.  AVERAGE  RADAR  CROSS  SECTION 

Important  developments  in  understanding  the  nature  of  sea  clutter  have  been  made 
by  studying  the  interaction  of  electromagnetic  waves  with  the  surface  of  the  ocean.  These 
studies  [2-61  ire  based  on  Rice’s  work  [1]  on  slightly  rough  surfaces  (surfaces  whose 
height  varia  ms  are  small  in  relationship  to  the  incident  wavelength  and  whose  slopes  are 
«1).  Us...g  perturbation  techniques,  they  found  that  the  reflected  energy  was  directly 
proportional  to  the  energy  density  spectrum  of  the  surface  height  variation  evaluated  at 
the  Bragg  scattering  condition.  Specifically,  the  radar  cross  section  for  direct  polarization 
is 

°HH  = sin40  cchjjW( 2k  cos 0)  (1) 

ovv  = 47rfe4  sin40  aVy\V(2k  cos 0) , (2) 

where  k = 2? r/A.  is  the  wavenumber  of  the  incident  wave,  0 is  the  grazing  angle,  W(K)  is 
the  energy  tensity  spectrum,  and  K-  2k  cos  0 is  the  Bragg  resonant  condition.  (Valen- 
zuela [4]  fwnd  the  cross-polarization  cross  section  by  considering  the  second-order  per- 
turbation terms.)  The  a terms  are 


and 


<*hh  = 


(e  -1) 


2 


[sin0  + (e  -cos20)1/2]2 


(Xyy  = 


(e  - l)[(e  - 1)  cos20  + ej 
[e  sin0  + (e  - cos2#)1/2]2 


(3) 


(4) 


where  e is  the  complex  dielectric  constant  of  the  ocean.  Now,  if  the  wave  spectrum  S(co) 
is  available,  W(K)  can  be  calculated  since 


iinV(X)  dK  = S(<o)  doo  (5) 

where  co2  = gk.  However,  a more  useful  approach  is  to  use  a result  of  Phillips  [11]  who 
investigated  the  growth  of  water  waves  and  concluded  that  there  exists  an  upper  bound 
for  the  height  of  gravity  and  capillary  waves.  Using  a dimensional  argument,  he  concluded 
that  in  the  range  where  waveheight  is  limited,  the  energy  spectrum  was  of  the  form 


W(tf)  = BIT4  . (6) 

There  is  uncertainty  about  the  value  of  B;  however,  Phillips  gives  B « 6 X 10" 3 for  gravity 
waves  and  B « 1.5  X 10-2  for  capillary  waves.  Substituting  Eq.  (6)  into  Eqs.  (1)  and  (2), 
one  obtains  the  following  limiting  values  for  the  cross  sections: 
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Ojjf-i  — 1 • £j7t  X 10  tan40 


(7) 


and 


avv  = 1.57T  X 10"3ayy  tan40  , 


(8) 


where  the  a’s  are  defined  by  Eqs.  (3)  and  (4). 


Guinard  and  Daley  [7]  have  compared  sea  clutter  data  with  the  theoretical  model. 
Their  data  were  taken  with  the  four-frequency  radar  (4FR)  system,  which  is  an  airborne, 
coherent,  pulsed  radar  capable  of  transmitting  a sequence  of  four  frequencies,  X band 
(8910  MHz),  C band  (4455  MHz),  L band  (1228  MHz),  and  P band  (428  MHz),  alterna- 
tively on  horizontal  and  vertical  polarizations.  The  basic  radar  parameters  of  the  4FR 
system  are  given  in  Table  1.  A basic  description  of  the  4FR  system  can  be  found  in 
Ref.  7 and  a detailed  description  in  Guinard  [12]. 


The  data  used  in  Guinard  and  Daley’s  study  [7]  came  from  two  major  experiments. 
In  July  1965,  a measurement  program  was  conducted  off  Puerto  Rico.  Radar  echoes 
were  recorded  for  all  frequencies,  both  polarizations,  and  a variety  of  grazing  angles  and 
sea  states.  Sea  conditions  were  measured  by  a team  from  the  Applied  Physics  Labora- 
tory, Johns  Hopkins  University.  During  the  measurements  the  sea  state  varied  between 
0 and  2 and  the  maximum  wind  was  20  knots  with  9-ft  seas.  The  second  measurement 
program  was  conducted  in  February  1969  in  the  North  Atlantic  to  study  rougher  sea 
conditions.  Ground  truth  was  provided  by  ocean  stations  India  (59°N,  19°W)  and  Juliet 
(52.5°N,  20°W).  Sea  states  varied  between  0 and  8,  a maximum  wind  of  45  knots  was 
recorded,  and  a significant  wave  height  of  26  ft  was  observed.  The  L-band  data  are  shown 
in  Figs.  2 and  3 for  vertical  and  horizontal  polarizations,  respectively.  Each  data  point 
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Fig.  2— Variation  of  RCS  with  grazing  angle;  W(K) 
e = 73-85i  (7] 
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DEPRESSION  ANGLE  (DEGREES) 

Fig.  3— Variation  of  RCS  with  grazing  angle-  W(K)  = 6 X 10“3K“4, 
e = 73-85i  [7] 

represents  the  median  value  of  cross  section  obtained  by  processing  approximately  30  s of 
data.  The  pulsewidth  selected  for  all  points  was  0.25  jus.  On  July  15,  19,  and  21,  the 
seas  were  3 to  5 ft  high,  on  July  20  they  were  2 to  3 ft  high,  and  on  July  27  and  28 
they  were  less  than  1 ft  high.  In  Fig.  2 (vertical  polarization),  one  can  see  that  the  theo- 
retical bound  of  Eq.  (8)  provides  a realistic  upper  bound  for  sea  clutter  and  exhibits  the 
correct  variation  with  grazing  angle. 

However,  the  situation  is  more  complex  when  the  data  are  compared  to  the  theoret- 
ical limit  in  Fig.  3 (horizontal  polarization).  V hile  there  is  good  agreement  for  large 
grazing  angles,  there  is  a wide  decrepancy  between  model  and  data  for  small  grazing 
angles.  To  understand  the  source  of  this  discrepancy,  recall  that  in  the  ocean  the  small 
waves  (i.e.,  Bragg  scatterers)  are  sitting  on  th°  much  larger  gravity  waves  and  swells.  In 
the  composite-surface  scattering  model  (wl,..,n  defines  the  role  of  the  large  and  small 
waves),  Wright  [3]  noted  that  the  only  significant  effect  of  the  large  waves  is  that  they 
change  the  angle  between  the  reflecting  surface  (ocean)  and  incident  radiation.  This  an- 
gular change  can  be  resolved  into  two  components:  one  a change  in  the  apparent  grazing 
angle;  the  other,  a rotation  normal  to  the  plane  of  incident  radiation.  Obviously,  the  first 
component  increases  cross  sections,  and  Valenzuela  [5]  has  shown  that  considerable  varia- 
tion in  cross  section  can  be  caused  by  the  latter  component.  Guinard  and  Daley  [7]  con- 
clude that  to  bound  the  cross  section  observed  with  horizontal  polarization,  it  is  necessary 
to  use  the  vertical  bound  of  Eq.  (8)  for  vertical  and  horizontal  polarization. 

The  purpose  of  this  section  was  to  introduce  the  slightly  rough-surface  and  composite- 
surface  scattering  models.  Those  with  a further  interest  in  the  subject  should  consult 
Guinard  and  Daley  [7] , which  provides  an  overview  of  the  theory  and  references  all  the 
important  work. 
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III.  DENSITY  FUNCTION  OF  SEA  CLUTTER 

In  this  section,  the  results  of  fitting  the  clutter  distribution  with  the  log-normal  and 
contaminated-normal  distributions  are  reviewed.  Then,  from  the  composite-surface  scat- 
tering model,  a conditional  density  p(x|ao)  is  introduced,  is  shown  to  be  Ricean,  and  is 
used  to  explain  the  non-Rayleigh  density  p{x)  of  sea  clutter.  Finally,  the  variation  of  the 
clutter  densities  p{x)  and  p(jc|ao)  with  parameters  such  as  radar  frequency,  pulsewidth, 
polarization,  and  wind  direction  is  found  by  analysis  of  variance  techniques. 


Fitting  of  Clutter  Data 
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Since  1951  [8] , it  has  been  known  that  the  density  function  of  sea  clutter  was  not 
Rayleigh  if  the  radar  pulsewidth  was  small.  While  some  research  was  done  by  laboratories 
and  industry  on  the  non-Rayleigh  density  (for  example,  Macdonald  [13]  and  Ballard  [14]), 
nothing  appeared  in  journals  until  Trunk  and  George  [10].  They  considered  fitting  the 
log-normal  and  contaminated-normal  distributions  to  the  distribution  of  envelope-detected 
sea  clutter.  (Through  the  remainder  of  this  report,  unless  otherwise  stated,  “distribution 
of  sea  clutter”  will  mean  the  distribution  of  the  return  of  envelope-detected  sea  clutter.) 

By  log-normal  clutter,  it  is  meant  that  the  envelope-detected  sample  x from  the  sea  re- 
turn has  the  density  function 


p(x)  = — exp 

( 2iro2x 2)1^2 


-2[in(xlxm)]‘ 


where  xm  is  the  median  value  of  x and  o is  the  standard  deviation  of  [fin  x] 2.  By 
“contaminated-normal  clutter”,  it  is  meant  that  the  quadrature  components  y have  a 
contaminated-normal  density 


(1  “ 7)  ( ~y^\  7 f ~y^ 

P(y)  ~ (2tt a2)1/2  exp\2a2/+  (2itK2o2)  eXP\2K2o2 


where  7 is  the  contamination  fraction  and  K is  the  ratio  of  the  standard  deviations  of  the 
two  Gaussian  densities.  In  a straightforward  calculation,  it  can  be  shown  that  the  density 
of  the  envelope-detected  sample  is 


p(x)  = (I-7)2 


27(1  - y)x 

+m^~eKp 


+ exp  

K2a2  \2K2a2 


-x2(K2  + 1)\  . (x2(K2  - 1) 
. 4 K2a2  i !°\  4K2o2 


Trunk  and  George’s  data  source  was  a frequency-agile,  high-resolution  radar  (FHR), 
which  is  an  airborne,  noncoherent,  pulsed  X-band  radar  capable  of  frequency  diversity  on 
a pulse-to-pulse  basis.  The  radar  operates  with  either  a long  (100-ns)  or  short  (20-ns)  pulse 
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with  vertical  or  horizontal  polariz-.tion.  The  beamwidth  is  0.5°  and  the  PRF  is  2,560  pps. 
Standard  operating  procedure  calls  for  the  aircraft  to  fly  at  180  knots  sampling  the  data  at 
a range  of  2 n.mi.  to  obtain  a grazing  angle  of  4.7°. 

Some  examples  of  a non-Rayleigh  sea  clutter  distribution  taken  by  NRL  in  1967  [15] 
are  plotted  in  Fig.  4 on  probability  paper  in  decibels,  with  the  median  of  the  density  arbi- 
trarily set  to  0 dB.  Each  curve  represents  a 2-min  sampling  interval.  As  one  can  see,  the 
sea  clutter  distribution  is  non-Rayleigh;  the  higher  the  sea  state,  the  more  non-Rayleigh  it 
becomes. 

The  clutter  models  were  fitted  to  the  data  by  a minimax  method.  First  the  param- 
eters, xm  of  the  log-normal  and  a of  the  contaminated-normal  densities,  were  used  to 
equate  the  medians  of  the  theoretical  distributions  with  the  median  of  the  actual  data. 

Then  the  remaining  parameters,  o of  the  log-noi.ual  and  y and  K of  the  contaminated- 
normal  densities,  were  used  to  minimize  the  maximum  difference  in  decibels  between  the 
theorized  curves  and  the  actual  data.  The  oest  fits  for  the  data  in  Fig.  4 are  shown  in 
Figs.  5 and  6.  Because  of  the  recording  methods  used,  data  recorded  before  1969  could 
not  be  thoroughly  processed. 

On  March  11,  1969,  the  FHR  radar  was  flown  200  mi  (320  km)  off  the  Virginia 
capes.  In  each  run,  the  aircraft  was  flown  in  a given  direction  with  respect  to  the  wind, 
i.e.,  upwind  (U),  downwind  (D),  or  crosswind  (C),  for  about  20  min.  After  2 min  either 
the  polarization  or  pulsewidth  was  varied  so  that  data  would  reflect  a variety  of  conditions. 


0 r.  4 S 8 10  12  14  16  18  20  22 


CLUTTER  CROSS  SECTION  IN  dB  TO  AN  ARBITRARY  SCALE 

Fig.  4— Experimental  sea  clutter  cross-section  data  taken  by  airborne  X-band  radar;  0.02-/is  pulse,  vertical 
polarization,  4.7°  grazing  angle.  Rayleigh  distribution  included  for  comparison  [10], 
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Data  on  sea  conditions  were  obtained  from  the  Fleet  Weather  Facility,  Naval  Reconnais- 
sance and  Technical  Support  Center,  Suitland,  Md.,  which  reported  25-to-31-knot  winds, 
8-ft  (2.5-m)  waves,  and  a 12-ft  (3.7-m)  swell. 


The  best  fitting  log-normal  and  contaminated-normal  distributions  for  some  short- 
pulse  (20-ns)  data  are  given  in  Table  2,  where  “best”  is  defined  by  the  minimax  solution. 
Columns  1 and  2 are  identifiers  which  distinguish  the  clutter  data  runs.  Column  3 is  the 
parameter  of  the  log-normal  which  yields  the  best  fit,  and  column  4 gives  the  maximum 
difference  (in  decibels)  between  the  data  and  the  best-fitting  distribution.  Columns  5-7 
contain  similar  information  about  the  contaminated-normal  distribution.  Both  distribu- 
tions provide  essentially  the  same  accuracy  of  approximation.  Also,  the  previously  re- 
ported observation  [9]  that  the  density  function  for  horizontal  polarization  has  a longer 
tail  than  that  for  vertical  polarization  seems  to  be  verified. 


Since  the  publication  of  Trunk  and  George  [10]  in  1970,  the  log-normal  density  has 
received  considerable  attention,  and  the  contaminated-normal  density  has  been  neglected. 
Before  we  leave  this  subject,  the  following  cautions  are  in  order: 


• Sea  clutter  is  not  log-normally  distributed.  (It  is  shown  in  the  next  section  to 
have  a spatially  varying  Ricean  density.) 


• While  the  log-normal  model  of  sea  clutter  can  yield  useful  results,  in  only  one  case 
does  it  yield  the  correct  probability  of  detection  in  clutter.  This  problem  is  discussed 
fully  in  Sec.  V. 


Table  2— Parameterization  of  Clutter  Data 


Clutter  Description 

Time  Identifier 


Log-Normal 

Maximum 

0 

Difference 

(dB) 

6.0 

1.6 

5.2 

1.6 

6.1 

1.3 

4.6 

1.6 

6.0 

1.0 

4.6 

0.4 

6.0 

1.0 

6.1 

0.8 

4.8 

0.8 

5.6 

1.8 

6.3 

0.8 

5.2 

1.2 

6.3 

0.8 

Contaminated-Normal 


K 

Maximum 

Difference 

(dB) 

0.025 

0.486 

0.034 

0.426 

0.051 


0.358 

0.104 
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B.  Density  of  Sea  Clutter 


If  the  density  of  sea  clutter  is  not  log-normally  distributed,  what  is  its  density?  Some 
insight  can  be  gained  by  examining  the  time  records  of  FHR  sea  clutter  data  shown  in 
Figs.  7 and  8.  Since  the  plane  is  traveling  at  180  knots,  each  record  corresponds  to  the 
the  return  from  an  approximate  2,000-ft  swath  of  ocean.  While  a structure  is  apparent  in 
both  sets  of  data,  it  is  much  more  obvious  in  Fig.  8.  The  basic  explanation,  provided  by 
Trunk  [16],  of  this  variation  with  time  is  that  the  width  (cr/2)  of  the  radar’s  illuminated 
patch  is  less  than  the  water  wavelength  of  the  sea  (over  60  m for  sea  state  5).  Conse- 
quently, the  density  of  the  envelope  return  x of  clutter  is  p[jc|a0(^g)] . That  is,  the  prob- 
ability of  any  value  x is  conditioned  on  the  average  backscatter  ao>  which  is  a function 
of  the  local  grazing  angle  gg  (Eqs.  (1)  and  (2))  at  the  range  cell  of  interest.  If  g is  the 
grazing  angle  for  a fiat  sea  and  s is  the  slope  of  the  large  wave  structure,  gg  = g + s.  Then, 
if  clutter  is  observed  over  a time  period  corresponding  to  several  water  wavelengths,  its 
density  function  can  be  written  as 
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p(x)  = fp[*la0(g  + s)]p[o0(s)]  da0(s) 
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Fig.  7— Sea  clutter:  short-pulse  data  taken  upwind  with  horizontal  polarization 


TIME  (s) 


« clutter:  short-pulse  data  taken  upwind  with  vertical  polarization 


where  p[a0(s)]  is  the  probability  density  of  a0,  which  is  a function  of  the  slope  of  the 
large  wave  structure  (swell  and/or  wind  waves).  Equivalently,  the  density  function  can  be 
written,  as 


P(x)  = fp[x\o0{g  + s)]p{s)ds. 


(13) 


where  p(s)  is  the  slope  density. 

For  simplicity  it  has  been  assumed  that  the  large  wave  structure  has  constant  slope 
over  the  illuminated  patch  area  A.  To  take  into  account  the  variation  of  slope  within  a 
range  cell,  let  S'(A)  be  the  variation  from  the  average  slope  s.  Then,  since  the  reflected 
power  from  elementary  Bragg  scattering  patches  in  a range  cell  add  noncoherently,  the 
average  cross  section  in  a cell  with  slope  s is  given  by 


0n(£  + s)  = T <f  a0[g  + s + S'(A)]  dA  , 


(14) 


NRL  REPORT  7986 


where  the  integral  is  a surface  integral  over  the  illuminated  patch.  Thus,  since  the  average 
slope  s is  a function  of  space,  the  non-Rayleigh  density  p(x)  is  seen  to  be  due  to  a spa- 
tially varying  density  p(tf|a0).  (In  reality  s is  a function  of  time  and  space.  However, 
since  the  data  are  taken  from  a plane,  the  major  change  is  spatial  and  for  our  purposes 
the  sea  may  be  considered  frozen.) 

Before  verifying  these  facts,  it  is  worth  noting  that  the  formulation  yields  a Rayleigh 
density  when  the  illuminated  patch  (specifically  the  length  along  the  wave  direction)  is 
large  with  respect  to  the  water  wavelength.  This  is  because  if  the  patch  encompasses 
many  waves,  the  density  of  the  average  slope  s is  a delta  function;  i.e.,  p(s ) = 8(s).  Con- 
sequently, p(x)  = p(*|o0);  an*4'  p(.tla0)  is  a Rayleigh  density  because  of  the  central 
limit  theorem. 


1.  Correlation  Properties 

For  large  illuminated  patches,  independent  samples  of  sea  clutter  can  be  obtained 
either  by  using  pulses  separated  by  about  10  ms  (the  decorrelation  time  usually  stated 
for  X-band  sea  clutter  return)  or  by  using  frequency  diversity.  However,  since  the  decor- 
relation is  due  to  phase  changes  between  capillary  waves,  for  a high-resolution  radar 
both  methods  should  yield  independent  samples  from  the  conditional  density  p(x|a0),  not 
from  the  density  p(x). 

To  verify  these  conclusions,  Trunk  [16]  calculated  the  correlation  functions  using 
12.8-s  intervals  (16,384  data  points)  of  data  taken  with  the  FHR  system.  As  can  be  seen 
from  Figs.  7 and  8,  this  interval  is  long  enough  to  represent  a sample  function  from  p(x). 
Decorrelation  times  are  given  in  Table  3.  (Decorrelation  time  being  defined  as  the  time 


Table  3— Correlation  Values  for  Data  Taken  by  FHR* 
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it  takes  the  correlation  function  to  fall  to  1/e.)  The  decorrelation  times  vary  from  12  to 
193  ms;  the  average  is  94  ms,  which  is  much  longer  than  the  quoted  time  of  10  ms. 

T?  ,'Se  particular  data  were  analyzed  because  the  pulse-to-pulse  frequency  diversity 
used  ' . each  12.8-s  interval  equaled  the  reciprocal  of  the  pulsewidth,  and  this  amount  of 
freq'  ency  difference  is  considered  sufficient  to  decorrelate  clutter  returns  [18] . The  cross- 
correlations for  the  two  frequencies  were  calculated  and  are  also  given  in  Table  3.  It  is 
seen  that  frequency  diversity  does  not  decorrelate  sea  clutter  returns  from  a high-resoluticn 
radar.  Next,  the  above  calculations  were  repeated  by  ividing  each  12.8-s  interval  into  64 
0.2-s  intervals  (short  enough  to  represent  a sample  function  from  p(xlo0))  and  averaging 
the  results  of  the  64  cases.  For  each  of  the  12  data  records,  the  decorrelation  time  was 
less  than  20  ms,  and  the  cross-correlation  was  less  than  0.1.  This  corresponds  favorably  to 
Pidgeon  [19],  who  reports  a maximum  correlation  of  0.2.  Thus,  time  separation  of  the 
samples  and  frequency  diversity  yield  independent  samples  from  p(xla0),  not  from  p(x). 


2.  Spatially  Varying  Ricean  Density 

What  is  the  density  function  for  the  conditional  density  p(ar|o0)?  First,  although  the 
illuminated  patch  is  rather  small  (10  ft  by  120  ft  (3  m by  36.6  m)  for  FHR  data),  there 
are  many  capillary  waves  in  it.  Consequently,  from  the  Central  Limit  Theorem  one  would 
expect  p(je|a0 ) to  be  Rayleigh  distributed.  However,  Trunk  [16]  has  shown  (the  analysis 
is  repeated  in  Appendix  A)  that  is  not.  Rather,  because  of  the  presence  of  dominant 
scatterers,  which  can  be  related  to  scattering  from  breaking  and  very  peak-crested  waves 
[20,21],  p(x|o0)  is  a Ricean  density. 

This  is  difficult  to  show  because  a0  changes  so  rapidly  (in  airborne  systems)  that 
there  are  too  few  samples  to  obtain  a good  estimate  of  p(x|o0)  before  o0  changes.  Fortu- 
nately, this  difficulty  can  be  overcome  by  analyzing  the  frequency-diverse  FHR  data  in  a 
special  way.  For  computation  ease,  a chi  density  will  be  used  in  place  of  the  Ricean 
density.  This  approximation  is  very  accurate  [22]  when  the  ratio  of  the  dominant  scat- 
terer  to  background  noise  is  in  the  neighborhood  of  0 dB. 

Let  and  yt  be  the  independent  samples  from  the  two  frequencies.  The  conditional 
densities  of  the  samples  are  then 


2xlf~1  exp  {-xf/2o2) 
p(x,\at)  r(k/2)(2a2)kl2 
and 


P(y,-I<7/) 


2yf'1  exp  (-yf/2o2) 

r(k/2)(2o2)kl2 


(15) 


(16) 


where  k is  the  number  of  degrees  of  freedom  (fe  = 2 is  the  Rayleigh  density)  and  Oj  is  a 
random  variable.  The  same  a,-  = a can  be  used  for  both  frequencies  because  the  time 
separation,  1/2560  s,  is  very  small  and  frequency  diversity  changes  a0  only  slightly. 
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As  mentioned  previously,  the  difficulty  in  showing  that  xt  and  y,  are  chi  (or  equiva- 
lently Ricean)  random  variables  is  that  o,  changes  too  rapidly  for  good  estimate  of  p(*,[a,) 
to  be  made.  Consequently,  a ratio  zi  = x(/y,-,  whose  density  is  independent  of  cq,  is 
formed.  The  density  of  z - z,  - xi/yi  is  found  by  forming  the  joint  density  of  x,-  and 
y;,  substituting  z = Xj/y,  and  yj  = yit  and  integrating  over  y This  yields* 


( \ = T(fe)  2 zk~1 

p(z)  r(fe/2)r(fe/2)  (1  + 22)fe 


(17) 


Using  the  FHR  data,  the  sample  distribution  of  independent  a,  (constructed  from  1,024 
samples  12  ms  apart,  i.e.,  every  16th  zt)  was  compared  to  Eq.  (17).  The  value  of  k was  ad- 
justed to  minimize  D,  which  is  the  maximum  difference  in  probability  between  the  sam- 
ple and  theorized  distributions.  The  fitting  results  are  shown  in  Table  4,  where  PD  is  the 
probability  that  the  maximum  difference  will  be  less  than  D when  the  theorized  distribu- 
tion, given  here  by  Eq.  (17),  is  the  true  distribution.  While  the  Kolmogorov-Smimov  one- 
sample  test  [23],  which  comparts  PD  to  a threshold  ft,  cannot  be  run,  since  k was  found 
by  minimizing  D (a  procedure  which  biases  the  test  in  favor  of  acceptance),  the  small 


Table  4— Fit  of  the  Chi  Density  to  the  FHR  Data* 


Identifier 

Optimum 

k 

D 

PD 

VUL 

2.7 

0.012 

0.003 

HUL 

2.4 

0.029 

0.663 

VDL 

2.4 

0.016 

0.068 

HDL 

2.6 

0.020 

0.237 

HDSt 

4.0 

0.035 

0.844 

VDS 

3.4 

0.021 

0.2S2 

VCS 

3.1 

0.029 

0.681 

HCS 

3.8 

0.019 

0.184 

VCL 

2.4 

0.021 

0.286 

HCL 

3.2 

0.017 

0.078 

HUS 

3.9 

0.014 

0.020 

VUS 

3.2 

0.017 

0.088 

♦From  reference  (16]. 

tThe  optimum  k for  this  case  is  greater  than  4.0.  The  search 
program  employed  search  in  the  interval  2 SftS  4.0  in  0.1 
steps. 
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♦The  reason  that  the  Ricean  density  was  not  directly  analyzed  is  that  the  density  for  z is  much  more  com- 
plicated than  Eq.  (17).  Specifically,  if  x and  y are  Ricean, 


P(2)  = 


2 z exp  (-R2/2) 


(1  + 22)2 

where  B is  the  ratio  of  the  dominant  scatterer  to  O. 
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values  of  D and  Pq  indicate  that  p(jc|o0)  is  at  least  nearly  a chi,  or  Ricean,  density.  For 
the  k’s  given  in  Table  4,  the  ratio  of  the  dominant  scatterer  to  background  noise  would 
vary  between  -3  and  +4  dB. 


3.  Dominant  Scatterers 

In  1974,  Lewis  and  Olin  [20]  measured  the  frequency  dependence  of  sea  return  with 
a system  that  transmitted  10-ns  simultaneous  pulses  at  8.6  and  9.2  GHz  from  the  same 
antenna  with  horizontal  polarization.  The  measurements  were  taken  at  the  Chesapeake 
Bay  Division  of  NRL  and  at  a very  short  range,  which  yielded  a range  cell  of  about  5 ft 
by  5 ft  (1.5  m by  1.5  m).  A typical  data  record  is  shown  in  Fig.  9;  the  waves  were  about 
4 ft  (1.2  m)  from  peak  to  trough,  with  whitecaps  forming.  Sea  return  was  found  to  have 
a large  dynamic  range,  with  the  largest  returns  coming  from  breaking  waves.  The  relation- 
ship between  sea  spikes  and  breaking  waves  was  recognized  and  confirmed  by  a foresight 
motion  picture  camera  mounted  on  the  radar  antenna.  The  camera  was  synchronized  to 
the  recorded  data,  viewing  a region  of  the  sea  containing  the  range  gate  responsible  for 
the  echoes.  It  should  be  noted  that  white  water  does  not  produce  sea  spikes.  Waves  must 
be  breaking,  a fact  which  suggests  that  spray  is  important. 
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The  fine  structure  in  the  sea  spike  was  studied  and  it  was  noted  that  the  return 
from  the  two  frequencies  appeared  to  be  uncorrelated.  Furthermore,  it  was  noted  that 
the  decorrelation  time  within  the  sea  spike  was  about  10  ms.  This  time  was  arrived  at 
by  noting  times  between  maximum  and  minimum  values.  Lewis  noted  that  these  facts 
could  be  explained  by  a fresnel  zoning  effect  produced  by  the  change  in  size  of  the 
whitecap  within  the  range  cell.  As  the  whitecap  changed  its  length  by  a fresnel  zone 
(A/4),  the  signal  amplitude  would  change  from  a maximum  to  a minimum  or  vice  versa. 
The  peak  return  was  approximated  by  assuming  that  a zone  had  an  effecti.e  height  of 
X/2  and  a width  w.  Then  the  scattering  cross  section  of  such  a zone  would  be 


-=  7 rw2  . 


The  cross  section  predicted  by  Eq.  (18)  was  compared  to  the  data  and  reasonable  agree- 
ment was  noted.  Research  is  still  being  done  in  this  area. 


M.  Long  [21]  has  also  noted  the  relationship  between  sea  spikes  and  large  wave 
structure.  He  reached  the  following  conclusions. 


1.  For  50%  of  the  sea  spikes,  a wave  breaks  (whitecap  forms)  simultaneously  or  a 
fraction  of  a second  thereafter. 


2.  About  40%  of  the  time,  a spike  was  called  when  a wave  structure  had  a very 
peaked  crest,  as  if  a whitecap  were  about  to  form  but  did  not. 


3.  No  breaking  waves  were  observed  in  the  absence  of  a sea  spike. 


Long  explains  the  sea  spikes  as  reflections  from  facets  (smooth  areas)  of  the  large  wave 
structure. 


While  the  mechanism  producing  the  sea  spikes  is  still  unknown,  there  is  no  uncer- 
tainty in  the  fact  that  sea  spikes  are  associated  with  breaking  waves  or  waves  that  almost 
break. 


Variation  of  Clutter  Densities 


The  scanning  rate  of  a search  radar  is  typically  6 to  15  rev/min.  Thus,  during  the 
scan  time,  the  large-scale  sea  structure  in  a range  cell  changes  very  little.  Consequently, 
if  the  illuminated  patch  is  smaller  than  the  water  wavelength,  the  radar  return  from  sea 
clutter  will  come  from  the  conditional  density  p(x|a0)  rather  than  from  p(x).  On  the 
other  hand,  if  scan-to-scan  processing  is  performed,  it  is  very  likely  (depending  on  scan 
rates,  radar  geometry,  and  sea  conditions)  that  samples  are  being  obtained  from  p(x). 

Since  the  density  will  determine  the  behavior  of  any  detector  used,  the  variation  of  clutter 
densities  p{x)  and  p(x|o0)  with  such  parameters  as  frequency,  polarization,  pulsewidth, 
and  wind  direction  will  be  analyzed. 
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1.  Variation  of  p(x) 


To  find  the  variation  of  p(x)  with  frequency,  the  4FR  data  taken  on  February  10, 
1969,  in  the  North  Atlantic  was  analyzed.  Specifically,  the  data  analyzed  by  Trunk  [16] 
were  taken  with  a 0.5-ps  pulsewidth  at  L and  X bands;  the  antenna  beamwidths  were  5° 
and  5.5°,  respectively;  horizontal  polarization  was  used  for  transmission  and  reception; 
and  the  pulse-repetition  frequency  (PRF)  was  683  pps.  The  aircraft  flew  at  200  knots, 
the  range  was  2,000  yd  (1,829  m),  the  depression  angle  was  10°,  and  the  azimuthal  angle 
(the  angle  between  the  radar  beam  and  the  wind  direction)  was  varied  from  0°  to  45°  in 
15°  increments.  Ground  truth  was  provided  by  Ocean  Stations  India  and  Juliet,  which 
reported  winds  between  30  and  35  mi/h  (48  and  56  km/h)  a sea  of  13.1  ft  (4.0  m),  and 
a swell  of  18  ft  (5.5  m). 


The  log-normal  density  (Eq.  (9))  was  fitted  to  the  data  by  equating  xm  to  the  sample 
median  and  then  finding  the  o that  minimized  the  maximum  difference  in  decibels  be- 
tween the  log-normal  distribution  and  the  sample  distmution.  The  log-normal  density 
was  used  instead  of  the  contaminated-normal  one  because  it  is  easier  to  interpret:  It  has 
only  one  parameter  (excluding  xm)  and  a larger  value  of  o indicates  a longer  tail  asso- 
ciated with  the  density.  A 36-s  sample  of  data  was  used,  the  data  were  fitted  between 
the  50  and  99.95  percentiles,  and  the  results  of  16  cases  appear  in  Table  5.  The  fitting 
errors  ranged  from  0.44  to  1.41  dB;  the  average  was  0.78  dB. 


From  Table  5 it  appears  that  the  clutter  distribution  is  a function  of  the  radar  fre- 
quency but  not  of  the  azimuthal  squint  angle.  However,  to  test  whether  these  factors 
cause  a significant  change  in  the  distribution,  a statistical  procedure  called  analysis  of 
variance  was  used.  This  procedure  decides  whether  differences  in  experimental  results  (in 
our  case  different  values  for  o)  are  true  differences  or  just  experimental  (sampling)  errors. 
The  results  of  the  analysis  of  variance  will  appear  in  this  section  and  further  details  are 
given  in  Appendix  B.  The  analysis  of  variance  procedure  was  applied  to  the  o’s  in  Table 
5 and  the  results  are  summarized  in  Table  6.  The  quantity  Sp  estimates  the  sum  of  two 
effects,  the  sampling  error  and  the  effect  of  changing  frequency;  the  quantity  S ^ also 
estimates  the  sum  of  two  effects,  the  sampling  error  and  the  effect  of  changing  azimuthal 
angle;  and  the  quantity  3'f  is  an  independent  estimate  of  the  sampling  error.  The  quantity 
c>p  is  compared  to  S|.  Since  their  ratio  is  largo,*  a frequency  effect  is  present.  Next,  S| 
is  compared  to  Sf . Since  their  ratio  is  small,  an  azimuthal  effect  is  not  present.  That  is, 
changing  the  azimuthal  angle  from  0°  to  45°  does  not  change  the  clutter  distribution 
significantly. 


P-band  data  (428  MHz)  were  fitted  to  the  log-normal  density  in  an  identical  manner: 
the  results  appear  in  Table  7.  The  P-band  distributions  have  a much  larger  spread  (higher 
o)  than  the  X-band  data  but  a slightly  smaller  spread  than  the  L-band  data.  Unfortu- 
nately, no  definite  comparison  can  be  made  between  these  bands  since  the  P-band  antenna 
had  a 12°  beamwidth  resulting  in  an  illuminated  patch  area  140  percent  larger  than  those 
of  L and  X bands.  However,  since  a larger  patch  tends  to  make  a distribution  have  a 


o o 

* Under  the  hypothesis  that  there  is  no  frequency  effect,  Sp  and  Sf  are  independent  estimates  of  the  sam- 


pl!ng  error,  and  their  ratio  has  an  F-distribution.  The  threshold  value  for  any  significance,  which  for  a 
significance  of  0.1  is  3.2,  can  be  found  in  either  Duncan  [24]  or  Fischer  [25],  which  also  provide  further 
information  about  analysis  of  variance.  All  analysis  of  variance  in  this  section  were  conducted  at  a 
significance  level  of  0.1. 
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Table  5 —a  in  Decibels  of  the  Fitted  Log-Normal  Distribution 


4FR  Radar 

Azimuthal  Squint  Angle 

Parameters 

0° 

j 15° 

30° 

45° 

Frequency 
L band 
X band 


♦From  reference  [16]. 


Table  6— Result?  of  Analysis  of  Variance:  The  Effect  of  Various 
Parameters  on  the  4FR  Clutter  Distributions* 


Effect 

Mean  Square 

Variance  Ratio 

Conclusion 

Frequency 
Azimuthal  angle 
Sampling  error 

= 5.640 

Sj  = 0.062 
Sf  = 0.101 

= 55.84 

Sj/Sf  = 0.61 

Frequency  effect  is  present 
Azimuthal  effect  not  present 

♦From  reference  [16]. 


Table  7— a in  Decibels  of  the  Fitted  Log-Normal 
Distribution  for  P Band* 

Azimuthal  Squint  Angle 


5.8  5.7  5.1  5.8  5.1  5.8  6.0  6.2 

♦From  reference  [16]. 
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smaller  spread,  P-band  distributions  have  a larger  spread  than  X-band  distributions,  but 
the  relationship  between  P and  L bands  cannot  be  determined.  Valenzuela  and  Laing  [26] 
analyzed  the  4FR  data  and  concluded  that  L-band  data  were  more  Rayleigh  then  X-band 
data;  based  on  results  of  K-S  tests.  However,  if  one  studies  the  results  of  the  individual 
teste,  one  sees  that  the  maximum  difference  in  probability  (D)  between  L-band  data  and 
the  Rayleigh  distribution  is  larger  than  the  difference  between  X-band  data  and  the  Ray- 
leigh distribution.  The  reason  Valenzuela  and  Laing  concluded  that  L-band  data  were 
more  Rayleigh  than  X-band  data  is  probably  because  Pq  is  smaller  for  the  L-band  than 
for  the  X-band.  However,  since  they  concluded  that  neither  data  are  Rayleigh  distributed, 
the  only  reason  PD  is  smaller  for  L-band  data  is  that  there  are  fewer  independent  samples 
for  L-band  data  than  for  X-band  data.  The  result  only  indicates  that  the  K-S  test  is  not 
a powerful  test;  i.e.,  if  data  do  not  come  from  the  theorized  distribution,  many  independ- 
p.it  samples  are  needed  to  reject  the  hypothesis  that  the  data  did  come  from  the  theorized 
distribution. 

The  FHR  data  are  analyzed  again  to  find  the  variation  of  p(x)  with  polarization, 
pulsewidth,  and  wind  direction.  First,  the  data  were  fitted  to  the  log-normal  distribution 
(results  appear  in  Table  8);  the  fitting  errors  ranged  from  0.4  to  1.4  dB,  and  the  average 
error  was  0.8  dB.  An  analysis  of  variance  was  conducted;  the  results  are  summarized  in 
Table  9.  The  most  significant  parameter  is  pola  ization:  the  clutter  distribution  for  hori- 
zontally polarized  data  has  a much  longer  tail  nan  that  for  vertically  polarized  data.  The 
next  most  significant  parameter  is  orientation  with  respect  to  the  wind:  upwind  or  down- 
wind measurements  have  a longer  tail  to  their  distributions  than  crosswind  measurements. 
While  the  analysis  of  variance  (Table  9)  shows  no  significant  difference  between  the  up- 
wind and  downwind  measurement,  an  analysis  of  variance  (Table  10)  of  the  short-pulse 
data  in  Table  2 definitely  shows  a difference  between  upwind  and  downwind  measure- 
ments. (Clutter  time  1145  was  eliminated  to  make  an  equal  number  of  cases  in  each 
category.)  Finally,  while  the  20-ns  data  have  a longer  tail  than  the  100-ns  data,  the 


Table  8— Log-Normal  Fit  of  Data  Taken  by  the  FHR* 


Date 

Time 

. 

Identifier 

a 

(dB) 

Maximum 

Difference 

(dB) 

3/11 

1344 

VUL 

4.8 

0.45 

3/11 

1346 

HUL 

6.0 

0.73 

3/11 

1411 

VDL 

4.7 

0.50 

3/11 

1413 

HDL 

7.3 

0.81 

3/12 

1026 

HDS 

7.5 

0.83 

3/12 

1028 

VDS 

5.5 

1.30 

3/12 

1049 

vcs 

4.4 

0.52 

3/12 

1051 

HCS 

6.1 

0.40 

3/12 

1125 

VCL 

4.4 

0.42 

3/12 

1127 

HCL 

5.6 

1.40 

3/12 

1154 

HUS 

5.9 

1.10 

3/12 

1156 

VUS 

5.8 

1.08 

♦From  reference  [16]. 
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Table  9— Results  of  Analysis  of  Variance:  The  Effect  of  Various 
Parameters  on  the  FHR  Clutter  Distributions* 


Source 

Mean 
Square  • 

Variance 

Ratio 

Conclusion 

Polarization 

6.75 

26.0 

Polarization  effect  is  present 

Pulse  length 

0.56 

2.2 

Pulse-length  effect  not  present 

Upwind  vs  downwind 

0.78 

3.0 

Effect  not  present 

Upwind  and  downwind 
vs  crosswind 

1.76 

6.7 

Effect  is  present 

Sampling  error 

0.26 

- 

♦From  reference  (16]. 


Table  10— Results  of  Analysis  of  Variance:  The  Effect  of  Various 
Parameters  on  the  FHR  Clutter  Distributions  in  Table  2 


Source 

Mean 

Square 

Variance 

Ratio 

Conclusion 

Polarization 

Upwind  vs  downwind 

Upwind  and  downwind 
vs  crosswind 

Sampling  error 

141.0 

4.9 

4;9 

Effect  is  present 
Effect  is  present 
Effect  is  present 

analysis  of  variance  (Table  9)  indicates  no  significant  difference.  Since  it  is  known  that 
the  clutter  distribution  is  Rayleigh  for  large  pulsewidth,  there  seems  to  be  a threshold 
effect.  That  is,  small  pulse  measurements  have  a larger  clutter  spread  than  large  pulse 
measurements;  however,  once  the  pulsewidth  is  smaller  than  the  water  wavelength,  changes 
are  no  longer  significant. 


2.  Variation  of  p(x|o0) 

To  find  the  effect  of  various  parameters  on  the  conditional  density  p(«|o0),  an  anal- 
ysis of  variance  was  run  on  the  data  in  Table  4 (optimum  k are  used  as  data  points);  the 
results  are  summarized  in  Table  11.  The  most  significant  parameter  is  pulsewidth:  the 
larger  the  pulsewidth,  the  more  Rayleigh  (as  oppose  to  Ricean)  the  density  of  clutter 
baekscatter.  This  can  be  explained  by  recalling  that  the  dominant  scatterers  are  breaking 
waves.  If  the  patch  is  small,  there  may  be  only  one  breaking  wave  in  it,  and  the  density 
will  be  Ricean.  However,  if  the  pulsewidth  is  increased,  several  breaking  waves  may  be 
present.  Since  return  from  these  waves  will  add  noncoherently,  the  density  will  tend 
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Table  11— Analysis  of  Variance  Results:  The  Effect  of  Various 
Parameters  on  the  Conditional  Clutter  Density* 


Effect 

Mean 

Square 

Variance 

Ratio 

Polarization 

0.606 

- 

7.88 

Pulse  length 

2.706 

35.21 

Wind 

0.006 

0.07 

Sampling  error 

0.076 

- 

.'-A  ••  V'-V 

L1/  .*=*.•> 


r • 


*From  reference  (16]. 
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toward  a Raleigh  density  as  the  pulsewidth  and  thus  the  number  of  breaking  waves  in- 
creases. The  other  significant  parameter  is  polarization,  vertical  polarization  giving  rise  to 
a more  Rayleigh-like  density  than  horizontal  polarization.  This  could  be  explained  by  the 
simple  fact  that  the  clutter  return  is  higher  for  vertical  polarization.  Thus,  if  the  dominant 
scatterers  have  the  same  cross  section  for  both  polarizations,  the  ratio  of  the  dominant  scat- 
terer  to  background  clutter  will  be  smaller,  and  hence  the  value  of  k will  be  nearer  to  2. 
Finally,  wind  direction  has  no  effect  on  the  conditional  density. 


D.  Empirical  Density  of  p{a0) 

In  1972,  Trunk  [16]  showed  that  the  non-Rayleigh  density  p(x)  of  envelope-detected 
sea  return  could  be  expressed  as 


p(x)  = J p(xlo0)p(o0)  do0 


While  Trunk  in*  astigated  the  densities  p{x)  and  p(x|oo)  and  concluded  that  the  form  of 
p(a0)  has  an  important  influence  on  p{x),  he  did  not  find  the  density  for  p(a0).  This 
situation  was  rectified  in  1974  when  Owens  [27]  used  a nonparametric  estimator  for  the 
density  of  10  log  o0,  which  is  equivalent  to  finding  p(o0).  Owens’  analysis  is  summarized 
in  the  following  sections. 

1.  Nonparametric  Estimation  of  a Probability  Density 

If  , . . . , aN  are  independent  random  variables  with  a common  density  f,  a kernel 
estimator  of  f at  a point  y has  the  form 


1 A 

f(y)  = ]v  I ] K(y'  2i> 


where  K is  a known  kernel.  A common  form  of  K is 
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Now,  consider  the  following  procedure.  For  each  i,  {i  = 1,  N],  select  a sample  ni  from  f, 
and  then  select  k samples,  denoted  xn,  . . . , xik  from  the  conditional  density  p(x|o,). 
Since 

E{Xjj}  =•  bOi  for  = 1,  ....  k,  (28) 


it  follows  that  an  estimate  of  bOf  is 


b5i  = \ E XU  ■ 

7 = 1 


Takx  ig  the  log  yields 


ir/NlCo.* 


p,-  = 10  log  bbi  = 10  log  dj  + 10  log  b . (30) 

That  is,  q{-  is  an  estimate  of  10  log  ot,  shifted  by  a constant  that  is  independent  of  i.  Thus, 
a i defined  by  Eqs.  (29)  and  (30)  can  be  used  in  Eq.  (23)  to  estimate  f(y). 


2.  Data  Analysis 

The  data  analyzed  by  Owens  were  the  small-pulse  (20-ns)  FHR  data  that  appear  in 
Table  8.  The  nonparametric  estimation  procedure  is  applied  to  these  data  in  the  follow- 
ing manner.  First  the  decorrelation  times  for  selected  short-pulse  data  are  extracted  from 
Table  3 (presented  in  column  3 of  Table  12).  Since  10  ms  is  the  decorrelation  time  of 
sea  clutter  when  o0  remains  constant,  the  decorrelation  times  in  Table  12  are  good  esti- 
mates of  the  times  required  for  o0  for  the  illuminated  patch  to  decorrelate.  Based  on 
these  times,  Owens  selected  a “sampling  interval”  (approximately  one-eighth  of  the  de- 
correlation  time),  in  which  it  is  assumed  that  a remains  constant.  For  each  data  record, 
every  eighth  sample  was  selected  to  form  a new  record,  with  samples  3 ms  apart  and 
adjacent  samples  recorded  at  different  frequencies.  The  total  number  of  samples  (this 
corresponds  to  k in  Eq.  (29))  in  a given  sampling  interval  is  given  in  column  5 and  the 
total  number  of  intervals  (corresponding  to  Not)  used  in  the  analysis  is  given  in  column  6. 
Then,  using  the  data  sets  indicated  in  Table  12,  Owens  computed  the  densities  as  follows: 


Table  12— Decorrelation  Times  and  Sampling  Intervals 


Identifier 

Decorrelation 

Time 

(ms) 

Sampling 

Interval 

(ms) 

Number  of 
Samples 
in  Interval 

Number  of 
Intervals 

HDS 

131 

39 

14 

208 

VDS 

168 

51 

18 

256 

VCS 

31 

9 

4 

192 

HCS 

112 

33 

12 

256 

HUS 

106 

33 

12 

256 

HUS 

193 

63 

22 

256 
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1.  Use  Eq.  (29),  to  obtain  average  the  returns  in  each  sampling  interval. 

2.  Compute  the  estimates  o(-  of  o0  (in  decibels)  for  each  interval,  using  Eq.  (30). 


3.  Estimate  the  probability  density  of  a0  (in  decibels),  using  Eq.  (23). 

The  density  of  each  of  the  six  cases  is  given  in  Ref.  27,  and  the  density  functions  for  the 
downwind  cases  are  shown  in  Figs.  10  and  11. 

The  density  functions  for  the  downwind  and  upwind  cases  are  quite  irregular.  This 
may  be  the  result  of  shadov/ing  of  the  patch  by  large  waves.  Furthermore,  densities  for 
vertical  polarization  are  more  peaked  than  those  for  horizontal  polarization.  This  sup- 
ports the  conclusion  in  (Ref.  16)  that  the  density  p(x)  has  a longer  tail  for  horizontal 
polarization  than  for  vertical  polarization.  Further  information  about  p(o0)  can  be  found 
in  Sec.  IV.B. 
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Fig.  10— Estimated  probability  density  of  normalized  RCS 
(horizontal  polarization,  downwind) 
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Fig.  11— Estimated  probability  density  of  normalized  RCS 
(vertical  polarization,  downwind) 


SEA  SURFACE  SIMULATION 


Trunk  [16]  and  Owens  [27,29]  have  found  that  a simulation  of  the  sea  surface  can 
be  extremely  useful  in  certain  problems  associated  with  high-resolution  radars.  In  this 
section,  this  simulation  will  be  discussed,  and  two  examples  of  its  use  will  be  given. 


Sea  Surface 


All  investigators  [16,27,29]  have  constructed  realizations  of  the  sea  surface  by  using 
the  method  suggested  by  Neumann  and  Pierson  [30].  The  realization  is  obtained  by  con- 
sidering the  linear  solution  to  the  Lagrangian  equations  of  motion  for  the  sea.  The  long- 
crested  waves  are  given  by  the  parametric  equations  in  5, 


r /o2 

x(t)  = 5 - sin  I — (6  cos0  +y  sin0)  - co,f  , 
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y(t)  = y , 


and 


2(f)  = cos 


or 


8 


(5  cos0  + y sin0) 


“ «,*]  . 


(31) 


where  (x,  y,  z)  are  the  coordinates  of  the  sea  surface  and  0 is  the  direction  in  which  the 
wave  crest  is  moving  with  respect  to  the  x-axis.  The  realization  is  constructed  by  form- 
ing a linear  sum  of  elementary  solutions,  i.e., 


N 


X(t) 

y(t) 

z(t) 


5 - I]  ai si 
1=1 


sm 


\-j-  (5  cos 0,-  +y  sin 0,)  - co,-t  + yt 


= y 


N 

L 

i=i 


co 


18 


~ (5  cos0,-  + y sin0,-)  - co,  t + 7 ,- 


(32) 


where  7,-  are  independently  distributed  phases  between  o and  2tt.  Given  a wave  spectrum 
S(co),  Trunk  [16]  suggests  setting  the  N frequencies  co,-  by 


fw/  2 i - 

I S(o)dco=4w 

J 0 


i/‘ 


S(co)  dco 


(33) 


'0  - JO 

and  letting  a,-  be  a Gaussian  random  variable  whose  variance  is 


a2  (a,-)  = 


l)  = "1 


S(co)  dco 


(34) 


The  wave  spectrum  that  was  used  was  the  Kitaigorodskii  [31]  spectrum  for  a fully  de- 
veloped sea, 


dg2 


S( co)  = — exp  [-6(g/uco)4]  , 
co5 


(35) 


where  d = 0.0081,  b = 0.74,  g is  the  acceleration  of  gravity,  and  u is  the  windspeed. 


To  see  the  type  of  realizations  this  method  yields,  let  us  construct  two  realizations 
for  a 20-knot  wind.  First,  N = 100  values  of  co,-  were  chosen  by  Eq.  (33),  and  a2(a,-)  was 
calculated  using  Eq.  (34).  Next,  the  wave  directions  0,-  were  chosen  from  two  Gaussian 
densities:  the  first  with  a standard  deviation  of  0.2  rad  and  the  second  with  a standard 
deviation  of  0.5  rad.  The  surfaces  are  shown  in  Figs.  12  and  13.  Each  surface  is  about 
700  ft  (213  m)  long  and  300  ft  (91  m)  wide.  The  40  y cuts  are  7.5  ft  (2.3  m)  apart  and 
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Fig.  12— Simulation  of  700-ft  by  300-ft  (213.4*  by  91.4-m)  sea  surface,  with  a(0,)  = 0.2 


Fig.  13— Simulation  of  700-ft  by  300-ft  (213.4-  by  91.4-m)  sea  surface,  with  0(0,)  = 0.5  [16] 


the  5 values  (hence  approximate  x values)  are  1.4  ft  (0.4  m)  apart.  From  appearance,  Fig. 
13,  with  o(0,)  = 0.5,  is  the  more  realistic  sea  surface.  (Unfortunately,  while  Trunk  [16] 
used  o(0)  = 0.5  he  reported  a (9)  = 0.2.  Consequently,  Owens  [27,29]  later  used  o(0,)  = 
0.2.  It  is  the  author’s  opinion  that  this  fact  will  have  little  effect  on  upwind  and  down- 
wind results  but  could  have  a larger  effect  on  crosswind  results.)  While  there  have  been 
several  measurements  of  directional  spectra,  the  results  of  one  of  which  is  reported  in 
Ref.  30,  Phillips  [32]  states  that  no  simple,  realistic  model  for  the  direction  spectrum, 
corresponding  to  the  existing  one-dimensional  wave  spectrum,  exists. 


There  are  several  ways  in  which  the  sea  surface  realization  can  be  used.  In  Sec.  IV.B 
it  is  used  to  estimate  the  densities  for  p(x)  and  p(o0),  and  in  Sec.  IV. C it  is  used  to  cal- 
culate the  probability  of  detecting  a small  target  on  the  surface  of  the  ocean. 
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B.  Estimation  of  the  Densities  p(x)  and  p(o) 

From  Sec.  II,  the  cross  section  for  horizontal  polarization  at  a grazing  angle  0 is 

°HH(0)  = 47rfc4  sin40  ^ — — W(2k  cosd)  i 

[sin0  + (e  -cos20)1/2]2 

where 


e is  the  complex  dielectric  constant 
k = 2n/\ 

KW(K)  dK  = S(w)  du> 
oo2  = Kg 
K = 2k  cos0 

S( co)  is  the  wave  spectrum  of  the  sea  surface. 

Since  the  reflected  power  from  elementary  Bragg  scattering  patches  in  a range  cell  add 
noncoherently,  the  average  cross  section  in  an  illuminated  patch  area  A is  given  by 


JJ  y,  z)  dx  dy 


where 


(x,y,z)eA 


z is  the  height  of  the  sea  surface  at  the  point  (*,  y), 
a = a(*,  y,  z)  is  the  local  grazing  angle  at  the  point  (x,  y,  z), 

0(a)  is  the  reflected  power  given  by  Eq.  (36),  and 

£(*>  y,  z)  is  the  normalized  two-way  antenna  power  gain  at  the  point  (x,  y,  z). 

Using  Eqs.  (12),  (36),  and  (37),  Trunk  [16]  calculated  p(x)  in  the  following  manner: 
First,  the  sea  surface  shown  in  Fig.  13  was  constructed.  Then,  0,-  (the  average  radar  cross 
section  for  the  ith  sample)  is  calculated  by  approximating  the,  surface  integral  in  Eq.  (37) 
with  a double  summation  involving  210  points:  twice  the  azimuth  beamwidth  is  divided 
into  21  radials  (each  with  a separation  of  0.1  beamwdith)  and  the  range  is  divided  into 
10  equally  spaced  ranges.  For  each  of  the  210  unshadowed  points  the  local  grazing  angle 
a is  calculated,  o^^(a)  is  calculated  using  Eq.  (36),  and  0,-  is  calculated  by  Eq.  (37). 

Next,  o,  + 1 was  calculated  by  advancing  the  sea  surface  10  ft  (3.0  m)  and  repeating  the, 
calculations.  This  process  was  repeated  until  50  values  of  of  were  calculated.  Then,  p(x), 
given  by 

P(x)  =J  p{x\o)p(o)da  (12) 

was  approximated  by 


ill 
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P(X)  ~ 


x exp  [~x2/2of] 


where  the  ergodic  theorem  has  been  used  to  replace  the  average  over  p(a)  by  the  time 
series  in  o,  and  p(*|a)  has  been  assumed  to  be  a Rayleigh  density.  The  Rayleigh  density 
has  been  used  instead  of  the  Ricean  density  because  it  is  easier  to  manipulate  and  Trunk 
[16]  has  shown  that  it  will  not  affect  the  results.  That  is,  the  variation  of  p(x)  from  the 
Rayleigh  is  basically  due  to  the  density  p(o).  Several  distributions  were  calculated  [33] 
and  are  shown  in  Fig.  14.  These  calculated  distributions  indicate  the  effects  of  various 
parameters  whose  significances  have  previously  been  demonstrated  (Sec.  III.C):  (a)  the 
clutter  distribution  for  horizontal  polarization  HUS  has  a longer  tail  than  that  for  vertical 
polarization  VUS;  (b)  the  distribution  for  the  short  pulse  VUS  has  a longer  tail  than  that 
for  the  long  pulse  VUL;  and  (c)  the  distribution  for  the  upwind  case  HUS  has  a longer 
tail  than  for  the  crosswind  case  HCS.  (The  computer  model  does  not  distinguish  between 
upwind  and  downwind.)  Also,  the  20-knot  case  has  a longer  tail  than  the  15-knot  case. 
The  Rayleigh  curve  was  obtained  for  HUS  when  the  windspeed  was  2 knots;  i.e.,  when 
the  range  cell  is  greater  than  the  water  wavelength,  the  model  yields  a Rayleigh  density. 


Owens  [29]  calculated  p(a)  in  a similar  manner.  Using  o(0t)  = 0.2  (Fig.  12),  he  cal- 
culated 200  values  of  o,-.  Then  Eq.  (23)  was  used  to  calculate  the  density  for  10  logo0. 
The  results  are  shown  in  Figs.  15  and  16.  In  comparing  the  experimental  results  with  the 
simulation,  the  following  paragraph  of  Owens  [29]  is  quoted, 


“. . . one  should  remember  that  the  observations  of  log  a q Ere  obtained  through  different 
mechanisms  in  the  two  cases.  In  the  simulation  the  only  errors  associated  with  observations 
of  log  Oq  are  caused  by  the  inaccuracies  of  the  model  and  the  mathematical  calculations 
therein.  Whereas,  in  the  experiment,  00  is  observed  with  an  error  having  essentially  two 
components.  The  first  component  of  error  arises  from  the  fact  that  on  a given  sample  the 
only  observable  variable  is  the  return  x,  x being  related  to  a through  the  density  p(x|o0 ). 
This  component  of  error  is  further  complicated  ir  that  various  properties  of  p(xlao),  which 
is  not  known  explicitly,  depend  on  experimental  parameters.  To  reduce  this  component  of 
error,  a number  of  observed  values  of  x are  averaged  to  obtain  an  estimate  of  0,  the  average 
being  carried  out  under  the  assumption  that  o remains  constant  over  the  appropriate  time 
interval.  Of  course,  this  assumption  is  not  precisely  correct,  thus  introducing  a second  com- 
ponent of  error  in  the  estimate  of  O.” 


Owens  goes  on  to  say  that  “considering  the  complexity  of  the  mechanism  producing  sea 
clutter  at  low  grazing  angles  and  that  the  two  sets  of  results  are  not  based  on  the  same 
set  of  observables,  the  authors  feels  that  the  agreement  is  good.”  It  is  worth  noting  that 
Owens’  largest  disagreement  [29]  occurs  for  the  crosswind  cases.  However,  if  a(0,)  = 0.5 
is  used  instead  of  o(0,)  = 0.2,  the  spread  of  p(o)  will  increase  and  better  agreement  will 
result. 


Probability  of  Detecting  Small  Surface  Targets 


Some  small,  fast  ships  of  tomorrow’s  Navy,  such  as  surface-effect  ships,  will  probably 
use  collision-avoidance  systems  to  avoid  such  ocean  debris  as  oil  drums  and  logs.  Conse- 
quently, Owens  [29]  investigated  the  probability  of  detecting  targets  at  a given  height  h 
above  the  local  sea  surface.  It  is  important  to  the  debris-detection  problem  that  these 
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Fig.  16— Estimated  probability  densities  of  normalized  RCS 
(vertical  polarization,  downwind) 


targets  can  he  smaller  in  geometric  size  than  the  larger  ocean  waves,  and  so  the  line  of 
sight  from  antenna  to  target  may  be  obstructed  by  surface  waves. 

This  shadowing  problem  was  investigated  using  the  sea  surface  simulation.  At  any 
point  (xq,  yo)>  the  target’s  height  at  a time  t was  fi0(£)  - h'+  z0(£),  where  z${t)  is  given 
by  Eq.  (32)  During  periods  when  the  line  from  the  point  [x0,  yo,  h0(t)]  to  the  radar 
antenna  intersects  the  ocean  surface,  the  target  is  “shadowed”  and  cannot  be  detected. 
(This  ignores  refraction  effects,  which  are  presently  under  investigation.)  Various  cases 
were  computed  by  Owens  and  a typical  case  is  given  in  Fig.  17.  For  example,  if  the 
antenna  height  is  75  ft  (22.9  m),  target  height  h is  1 ft  (0.3  m),  and  the  sea  state  is  3 
(wind  about  15  knots),  the  probability  of  having  a clear  line  of  sight  at  any  moment  is 
0.5  at  a range  of  2 n.mi. 

Owens  also  calculated  the  probability  of  detecting  a target  by  calculating  the  clutter 
return  in  the  same  range  cell  as  the  target.  For  purposes  of  illustration,  consider  the  radar 
to  have  a beamwidth  of  1.5°,  a pulsewidth  of  20  ns,  and  a rotation  rate  of  60  rpm.  The 
radar  is  75  ft  (22.9  m)  above  the  surface,  aboard  a ship  traveling  at  80  knots. in  sea  state  3. 


V *- 


is® 


FiSeSs 


V £ »!?***? 

ppp. 

$rA>& 


mm 

C3 


mm 


CvvO-i 

(XT 


AAV 

mm 

iSP'-.’s: 


iSlf 

/%*  Vt ".-NV- 


1 %*  v v 

3 • * ■ * * * 


NRL  REPORT  7986 


hi 


bX 


£ 


t'-\ 


6 


r' 


i>  .*  • * - • 


Fig.  17— Probability  of  sighting  vs  range  for  an 
antenna  height  of  75  ft  (22.9  m)  and  a sea  state 
of  3 
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In  Fig.  18,  the  average  received  signal-plus-noise  and  average  clutter  level  are  graphed  as  a 
function  of  range. 


These  values  are  for  a particular  realization  of  the  sea  surface  and  both  curves  are 
referenced  to  the  ambient  noise  level.  First,  notice  the  lower  envelope  of  the  clutter  re- 
turn. This  represents  the  clutter  received  through  the  sidelobes;  when  the  clutter  level 
falls  to  this  envelope,  the  entire  illuminated  patch  is  obscured  by  large  wave  structure. 
The  upper  envelope  of  signal-plus-noise  indicates  the  target  return  when  the  target  is  not 
obscured.  When  the  curve  chops  away  from  the  envelope,  the  target  is  obscured  (zero 
effective  cross  section). 


Furthermore,  the  signal  and  clutter  returns  are  highly  correlated.  When  the  target  is 
not  obscured,  it  is  likely  to  be  on  the  front  side  or  top  of  a large  wave;  and  a large  clutter 
return  can  be  expected.  On  the  other  hand,  when  the  target  is  obscured,  the  illuminated 
patch  is  likely  on  the  back  side  of  a wave,  and  a relatively  small  clutter  return  will  result. 
Owens  assumed  a square-law  detector*  and  performed  scan-to-scan  integration.  Typical 
results  are  shown  in  Fig.  19.  As  will  be  seen  in  the  next  section,  substantial  improvement 
results  from  scan-to-scan  integration.  For  specific  details  of  this  work,  see  Ref.  31. 


V.  DETECTION  OF  TARGETS  IN  NON-RAYLEIGH  SEA  CLUTTER 


Since  1947,  the  classical  works  of  Marcum  [35]  and  Swerling  [36]  have  been  used 
to  calculate  the  probability  of  detecting  targets  in  sea  clutter.  Their  models,  based  on  a 


♦Since  a target  can  be  shadowed  for  a number  of  scans,  the  detector  incurs  a loss  similar  to  a collapsing 
loss:  A square-law  detector  is  used,  since  the  collapsing  loss  is  greater  for  a linear  detector  than  for  a 
square-law  detector  [34]. 
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non-Rayleigh),  in  order  for  the  target  to  remain  in  the  same  range  cell  it  must  either  be  almost 
stationary  or  some  kind  of  range  tracking  (during  detection)  or  collapsing  must  be  applied. 

Besides  the  problem  of  independence  of  the  samples,  there  is  the  question  of  the  ap- 
propriate density  for  the  clutter  return.  From  Fig.  20,  it  is  obvious  that  the  signal  required 
to  obtain  a particular  performance  level  is  very  dependent  on  the  distribution  assumed  for 
the  clutter— a log-normal  distribution  requires  more  signal  strength  than  a contaminated- 
normal  one,  which  requires  more  than  a Rayleigh  distribution.  Also,  the  differences  are 
larger  for  1 pulse  than  for  30  pulses.  Thus,  the  choice  is  extremely  critical  when  only  one 
or  a few  independent  samples  are  available.  It  is  this  author’s  opinion  that  for  small  false 
alarm  rates  (i.e.,  below  10-6),  the  log-normal  density  will  yield  pessimistic  results,  i.e.,  per- 
formance will  actually  be  better  than  predicted.  This  is  because  p(x ) is  not  really  log- 
normally  distributed  and  has  only  been  used  to  approximate  the  density  to  its  99.99  per- 
centile. As  shown  in  Sec.  III.B,  p(x)  eventually  falls  off  as  a Ricean  density;  i.e., 
exponentially.  Thus,  the  threshold  need  not  be  set  quite  as  high  as  the  log-normal  den- 
sity requires. 

Keeping  in  mind  the  previous  comments,  detection  results  are  presented  for  the  log- 
normal and  contaminated-normal  densities. 


A.  Log-Normal  Density 

For  the  log-normal  model  the  envelope-detected  sample  has  the  density  function 


p(x)  = 


( 2tro2x 2)1/2 


exp 


-2 [in  (x/xm)] 


(9) 


where  xm  is  the  median  value  of  x and  a is  the  standard  deviation  of  (in  a:)2.  Following 
Rice’s  procedure  [40] , we  can  show  [41]  that  the  density  function  of  a constant  signal  A 
in  log-normal  noise  is 


Fig.  20— Comparison  of  the  Rayleigh, 
coniaminaled-normal  (7  = 0.25,  K = 
2.25),  and  log-normal  (0  = 6 dB)  de- 
tection probabilities  for  N = 1 and  30 
pulses  and  Pfa  = 10“8.  For  this  com- 
parison, SjN  is  defined  in  terms  of  the 
median  value  for  all  the  probability 
distributions  (10], 
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t,A)=  f 

*Sf\ 


2n  x exp  [-2a-2  £n2(.x2  - 2xA  cosd  + A2)1^2]  dd 
7r(27ra2)1^2(jc2  - 2 xA  cos 6 + A 2) 


where  the  return  has  been  normalized  by  median,  so  that  xm  = 1. 

Various  detectors  will  now  be  evaluated  using  the  log-normal  model. 


1.  Mean  Detector 


The  straightforward  method  of  determining  the  probability  density  for  the  sum  of 
N independent  pulses  (i.e.,  the  mean  detector)  is  the  characteristic  function  method.  If 

0(co)  = E {exp(/co*)}=  I p(x.  A)  exp  (icox)  dx  (4( 

Jo 

is  the  characteristic  function,  the  probability  density  of  the  sum  of  N statistically  inde- 
pendent envelope-detected  pulses  is 

pN(x,  A)  = -^  f [<P(oj)]N  exp  (-icox)  doo  . (41 

The  probability  of  false  alarm  Pfa  for  a threshold  T is  given  by 

Pfa  = f Pn(x - 0)  dx  , (41 

Jrn 
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and  the  probability  of  detection  PD  is 


%-J' 


pN(x,  A)  dx  , 


With  the  fast.  Fourier  transform,  it  is  fairly  simple  to  calculate  Eqs.  (40)  and  (41). 
Then,  Pfa  and  Pq  are  calculated  by  straightforward  numerical  integration.  Threshold 
values  for  N = 1,  3,  10,  and  30  for  values  of  Pfa  ranging  from  10-2  to  10~8  are  shown 
in  Fig.  21  for  the  log-normal  density,  with  o = 6 dB.*  With  these  thresholds,  Eq.  (43)  is 
evaluated,  yielding  P^  vs  signal-to-noise  ratio  S/N  per  pulse,  where  S/N  is  the  signal  divided 
by  the  median  value  of  clutter.  This  is  not  the  usual  definition  of  S/N,  in  which  the  noise 
reference  is  the  rms  value  of  the  noise.  However,  it  is  an  appropriate  choice,  since  it  is 
the  median  value  of  clutter  that  is  usually  reported  [7].  The  PD  curves  for  a = 6 dB  are 
plotted  in  Fig.  22.  An  inspection  of  these  curves  shows  that  the  S/N  required  for 

♦Additional  threshold  values  and  detection  probabilities  for  a = 3 dB  and  0 = 9 dB  can  be  found  in  [41]. 

It  should  be  denoted  that  for  all  detectors  a is  not  known  a priori.  Consequently,  an  adaptive  thresholding 
technique  must  be  used. 
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detection  decreases  rapidly  as  N increases.  In  fact,  the  integration  gain  Ig  is  greater  than 
would  be  obtained  for  the  coherent  integration  of  a Gaussian  density,  i.e.,  Ig>  10  log  N. 

It  has  been  pointed  out  [37]  that  this  means  the  optimal  beamwidth  for  a high- 
resolution  scanning  radar  need  not  be  the  smallest.  For  instance,  consider  a case  in  which 
S/N  is  20  dB  per  pulse,  and  a three-pulse  scan-to-scan  integration  is  being  performed.  If 
the  clutter  has  a log-normal  distribution  with  a = 6 dB,  Pj)  will  be  0.1  for  a Pfa  of  10“6. 
If  the  beamwidth  and  scanning  rate  are  both  increased  by  a factor  of  10/3,  then  10  pulses 
can  be  integrated  in  the  same  time.  In  this  case,  while  S/N  falls  to  about  15  dB  per  pulse, 
PD  is  greater  than  0.99.  Thus,  better  performance  is  obtained  with  the  larger  beamwidth. 


2.  Binary  Integrator  and  Rank  Detector 

It  is  well  known  that  the  binary  integrator  (sometimes  called  a dual-threshold  detector 
or  Af-out-of-iV  integrator)  is  exactly  equivalent  to  a rank  detector.  Given  a set  of  N ordered 
samples,  Xi  < ...  < xk  < ...  < xN>  the  feth  ordered  rank  detector  involves  simply  com- 
paring xk  to  a threshold.  The  median  detector  is  the  special  case  where  k = N/2.  How- 
ever, since  the  implementation  and  evaluation  of  this  detector  is  simpler  in  the  form  of 
the  binary  integrator,  the  detector  will  be  discussed  in  this  form.  A simple  block  diagram 
of  the  binary  integrator  (or  ftth  rank)  is  shown  in  Fig.  23.  The  probability  p(A)  that  the 
returned  signal  exceeds  T is 


PW)  = f 


p{x , A)  dx 


where  p(x,  A)  is  given  by  Eq.  (39).  The  probability  that  more  than  m = N-  k of  the  N 
returns  exceed  T is  given  by  the  binomial  distribution 
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23— Implementation  of  binary  integrator  or, 
equivalently,  rank  detector  [10] 
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The  desired  probabilities  are  given  by  Pfa  = P(A  = 0)  and  Pq  = P(A).  In  1970  [10] , the 
median  detector  for  the  log-normal  distribution  with  a = 6 dB  was  evaluated.  The  thresh- 
olds are  presented  in  Fig.  24,  and  the  detection  results  are  presented  in  Fig.  25.  Com- 
paring the  mean  and  median  detectors  in  Fig.  26,  one  concludes  that  the  median  requires 
a smaller  S/N  than  the  mean  and  that  the  difference  increases  with  increasing  N:  for 
N = 3,  the  median  is  1.6  dB  better  than  the  mean;  and  for  N = 30,  the  median  is  2.8  dB 
better.  The  performance  difference  between  the  two  detectors  for  large  sample  sizes  can 
be  obtained  from  Pitman’s  asymptotic  relative  efficiency  (ARE)  [42] . This  criterion  is 
quite  appropriate  when  the  signals  are  very  weak.  Specifically,  it  is  the  ratio  of  the  num- 
ber of  samples  required  to  maintain  a Pfa  and  PD  for  the  first  detector  to  the  number  re- 
quired for  the  second  detector  as  the  S/N  approaches  zero: 


N^P^Po'S/N) 

ARE(d1,d2)  - ^ 


">o  %: 


where  N • is  the  minimum  number  of  observations  required  for  detector  d,-.  The  simplest 
way  of  calculating  the  ARE  is  by  employing  the  concept  of  efficacies,  which  was  also 
introduced  by  Pitman.  Using  a result  of  Noether  [43] , Trunk  [44]  calculated  an  ARE 
of  129.  (See  Appendix  C.)  Assuming  that  a noncoherent  integration  gain  of  10  log y/W 
is  appropriate  for  small  signals,  the  median  detector  is  approximately  10.5  dB  better  than 
the  mean  detector  for  very  small  signals. 
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Fig.  24— Normalized  bias  values  for  the  median  detector 
(log-normal  distribution  with  0 = 6 dB) 
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Since  the  log-normal  density  is  unsymmetrical  about  its  median  value,  there  is  no 
reason  to  expect  that  k = N/2  (the  median)  is  the  optimal  rank.  Rather,  because  of  the 
long  tails  of  the  log-normal,  one  would  expect  that  the  smaller  ranks  ( k < N/2)  would 
behave  better.  While  Trunk  [45]  showed  that  1 dB  could  be  gained  with  respect  to  the 
median  by  using  k = N/ 3 against  the  contaminated-normal  density,  Schleher  [38]  showed 
that  larger  gains  were  obtained  with  smaller  ranks  against  the  log-normal.  Schleher  [46] 
obtained  the  optimal  value  of  m for  different  values  of  N by  calculating  the  detection 
curves  for  each  possible  value  of  m and  choosing  the  best  value.  The  optimal  values  of  m 
for  N = 3,  10,  and  30  pulses,  for  the  log-normal  density  with  a = 6 dB,  are  given  in  Table 
13  along  with  the  optimal  values  for  the  Rayleigh  density,  which  were  determined  by 
Schwartz  [47] . Detection  curves  were  generated  by  Schleher  [46]  for  the  log-normal 
density  (a  = 6 dB)  with  N = 3,  10,  and  30  pulses.  The  detection  curves  for  the  optimal  m 
are  given  in  Figs.  27,  28,  and  29.  Comparing  these  results  with  the  median  detector  shows 
that  the  largest  differential  is  for  smaller  N.  For  N = 3,  the  median  is  about  4 dB  worse 
(for  PD  « 0.9  and  Pfa  = 10"6)  than  the  binary  integrator  (m  = 2);  for  N = 30,  the  median 
is  only  about  1 dB  worse. 
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Table  13— Optimal  Value  of 
Second  Threshold  m 

m (Log-Normal)  m (Rayleigh) 
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Fig.  27— Probability  of  detection  vs  S/N  for  the  log-normal  distribution 
(0  = 6 dB)  and  optimal  binary  integrator:  JV  = 3,  m = 2 
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However,  if  one  can  bound  the  performance  of  the  optima!  detector  and  find  a sub- 
optimal  detector  that  approaches  this  performance,  then  there  is  no  reason  to  specify  the 
optimal  detector. 

The  Chemoff  [53]  bound, 

Pfa  < exp  [p(s)  - sfi(s)] 

0 < s < 1 (50) 

PD  > 1 ~ exp  [n(s)  + (l-s)p(s)]  , 
where  the  semi-invariant  p(s)  is  given  by 

/oo 

[pN(x>A)]s[pN(x,A  = 0)]1~sdx,  (51) 

oo 

provides  a rather  coarse  bound  for  the  desired  probabilities.  Van  Trees  [54]  tightened 
the  bounds  by  finding  a multiplicative  factor  for  the  exponentials,  in  Ref.  50,  using  a 
central-limit-theorem  argument.  Schleher  [38,46]  obtained  a more  accurate  result  by 
using  an  Edgeworth  series  expansion  to  approximate  the  desired  probabilities.  Schleher 
[46]  shows  that  Pfa  is  given  by 


Pfa  = 


erfc  [a(S)]  exp  [a2(S)  + B(S)] 


exp  im]  {! erfc  [a(S)]  Is' v®] 3 exP 

[i-sms)]. 


where 


oc{S)  = SMSm1'2  (63) 

B(S)  = n(S)  - Sii(S) . (54) 

An  expression  for  1 - PD  can  be  found  [46]  by  substituting  1 - S for  S in  Eqs.  (52), 
(53),  and  (54).  Evaluation  of  Eq.  (52)  for  Pta  and  PD  requires  the  determination  of  the 
first,  second,  and  third  derivatives  of  (x(S).  Schleher  calculated  these  derivatives  by  eval- 
uating p(f ) at  many  points  in  the  interval  (0, 1),  fitting  a cubic  spline  function  to  the 
points,  and  differentiating  the  spline  function.  (The  AT-degree  spline  function  produces 
the  smoothest  curve  through  the  data  points  while  maintaining  continuity  of  the  first 
N - 1 derivatives  at  each  data  point.)  Receiver  operating  curves  were  generated  and  are 
given  in  Ref.  46.  Curves  comparing  the  mean,  median,  binary  integrator,  and  trimmed- 
mean  detectors  with  the  optimal  detector  are  given  in  Figs.  31,  32,  and  33,  for  N = 3, 

N=  10,  and  N = 30  pulses,  respectively.  From  these  curves  it  can  be  seen  that  the  binary 
integrator  and  trimmed-mean  detector  are  within  1 dB  of  the  optimal  detector.  Since  the 
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Fig.  33— Comparison  of  various  detectors  in  log-normal 
(a  = 6 dB)  interference 


trimmed-mean  detector  requires  a ranking  of  the  samples,  the  binary  integrator  should  be 
used  to  detect  nonfluctuating  targets  in  the  presence  of  log-normal  interference. 


5.  Fluctuating  Targets 


In  the  previous  sections,  PD  curves  were  generated  for  various  detectors.  It  has  been 
assumed  that  the  target  is  not  fluctuating  and  that  the  noise  samples  are  independent. 
However,  to  obtain  independent  samples  from  sea  clutter  (whose  density  is  approximated 
by  the  log-normal),  the  returned  samples  must  be  from  different  scans.  Since  the  time 
separation  is  probably  larger  than  1 s,  it  is  quite  likely  that  the  target  is  fluctuating.  Con- 
sequently, Trunk  [55]  calculated  the  performance  of  the  mean,  median,  and  trimmed- 
mean  detectors  against  the  log-normal  density,  with  Swerling  II  and  IV  fluctuations  as- 
sumed [36]. 


Since  the  threshold  value  for  any  Pfa  is  independent  of  the  target  model,  they  have 
been  already  calculated;  all  that  remains  is  to  calculate  the  PL  for  various  S/N.  Since  the 
interesting  range  of  PD  is  from  0.01  to  0.09,  Monte  Carlo  techniques,  which  require  con- 
siderably less  computer  time  than  the  characteristic-function  approach,  will  be  used.  It  is 
well  known  that  the  ith  sample  xt  of  an  envelope  detector  can  be  generated  by 


- (; * t 


+ Z 


>\l/2 


i) 


(55) 


where  y,-  and  2,-  are  in-ph2se  and  quadrature-phase  components, 
log-normal  clutter 


For  the  generation  of 
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exp  (a[-2  fin  (v^)]1^2  sin  (27ru,-2)}  sin  (2nvi3)  + A , 
exp  {cr[— 2 fin  (v/i)]1^2  cos  (27rt'i-2 )}  cos  (2iwi3) , 


(56) 


where  {v-j}  are  independent  and  uniformly  distributed  random  numbers  on  (0,  1)  and  A 
is  the  signal.  For  the  fluctuating  signal,  the  probability  densities  are 


p(A)  = M-  e~AZl°2 
a2 


(57) 


for  the  Swerling  II  case  and 


p(A)  = ~ e-2A2/a2 


(58) 


for  the  Swerling  IV  case.  The  generation  of  a variable  having  the  density  given  in  Eq.  (57) 
is  straightforward.  Integrating  Eq.  (57),  one  obtains 


P(u<A)  = 1 - e-^2/°5 


or 


e-A2/a2  = 1 - P(U<  A). 


(59) 


Now,  the  quantity  1 - P{u  < A)  is  uniformly  distributed  between  0 and  1 and  so  can  be 
replace^  by  a random  number  u.  If  Eq.  (59)  is  solved  for  A,  giving 


A = o[-fin(u)]1/2 


(60) 


then  A is  seen  to  have  a Swerling  II  distribution.  The  power  density  for  a Swerling  IV 
case  is 


P(Z)  = ~ e-2 Zlo  . 
a 2 


(61) 


If  7 - Xi  + jc2  and  if  p(x)  = (2/a)  exp  (~2x/a),  Z will  have  the  density  given  in  Eq.  (61). 
Consequently,  if 


A = 


-§  [fin  (ut)  + fin  (u2)] 


1/2 


(62) 


A will  have  a Swerling  IV  distribution.  With  Eqs.  (55),  (56),  (60),  and  (62),  a Monte 
Carlo  simulation  was  performed.  The  results  are  given  in  Ref.  55. 


In  Fig.  34,  the  fluctuating  results  are  compared  with  the  nonfluctuating  results.  This 
figure  illustrates  the  fact  that  for  a given  number  of  samples  N,  the  fewer  the  samples  in- 
volved in  the  detector  (i.e.,  n2  - nj  + 1)  the  larger  the  required  increase  in  S/N  to  maintain 
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the  desired  PD.  For  instance,  to  maintain  Pfa  = 10-6  and  PD  = 0.9  for  N = 3 and  a 
Swerling  II  target,  the  mean  detector  (which  uses  all  three  samples)  requires  a 4.6-dB  in- 
crease, whereas  the  median  detector  (which  was  only  one  sample)  requires  a 6.0-dB  in- 
crease. Similarly,  for  N = 10,  the  mean  requires  a smaller  S/N  increase  than  the  trimmed 


Since  the  binary  integrator  uses  only  ope  sample,  like  the  median,  this  author  be- 
lieves it  will  suffer  a “fluctuation  loss”  similar  to  that  of  the  median.  Thus,  for  fluc- 
tuating targets  and  log-normal  density,  the  trimmed-mean  detector  is  the  most  effective 
one,  the  binary  integrator  is  next,  and  the  median  is  third,  only  slightly  better  than  the 
mean.  Of  course,  for  large  samples,  the  behavior  of  the  detectors  for  fluctuating  targets 
will  approach  the  nonfluctuating  behavior.  Consequently,  for  large  sample  sizes,  the 
binary  integrator  and  trimmed-mean  detector  are  equivalent,  and  both  are  better  than  the 
median,  which  in  turn  is  better  than  the  mean. 


I A *1*  A «> 
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B.  Contaminated-Normal  Density 

For  the  contaminated-normal  model,  the  envelope-detected  sample  x has  the  density 
function 


p(*)  = (1-7) 


JL 


-x2\  + 72*2 


2a2/  K2o2 


27(1  - y)x 


-x2(K2  + 
, 4 K2a2 


m 


2 K2a2, 


x2(K2  - 1)^ 
, 4 K2a2  I 


where  7 is  the  contamination  fraction  and  K is  the  ratio  of  the  standard  deviations  of  the 
two  underlying  Gaussian  densities.  Following  Rice’s  procedure  [40].  we  can  show  [56] 
that  the  density  function  of  a constant  signal  A in  contaminated-normal  noise  is 


where 


2tt 

c,  A)  = f p{x,A,d)de  , 


p(x,  A , 0)  = — i (1  - 7)2  exp  [-(*2  - 2xA  cosff  +A2)/2a2] 
2ira2  { 


► . *.  *■  » .-  • « 
f * • i • 

* *'«  _ -■*  »L. 


+ — exp  [-(*2  - 2xA  cos 0 + A2)/2K2o2] 
K2 


+ — exp  [-(*2  sin20  - 2xAK2  cos0  + A2K2 
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+£2*2  cos20)/2£2a2] 


+ ^ exP  [-(*2  cos0  - 2*A  cos0  + A 2 

+£2x2  sin20)/2£2a2]l  . (64) 


Various  detectors  will  now  be  evaluated  using  the  contaminated-normal  model. 

1.  Mean  Detector 

The  threshold  T for  a desired  Pfa  and  the  pD  have  been  calculated  using  the  charac- 
teristic function  approach  given  by  Eqs.  (40)  through  (43).  The  thresholds  for  y = 0.25 
and  K = 2.25  can  be  found  in  Ref.  10,  and  the  detection  curves  for  N - 1,3,  10,  and  30 
pulses  are  given  in  Fig.  35,  where  S/N  is  still  the  signal  divided  by  the  median  clutter, 
which  for  this  case  is  1.41a.  Two  points  should  be  noted:  (a)  the  integration  gain  for 
the  contaminated-normal  density  is  greater  than  for  the  Rayleigh  density,  but  less  than  for 
the  log-normal  density;  (b)  comparing  Figs.  22  and  35  shows  that  the  contaminated- 
normal  model  requires  considerable  less  signal  strength  then  the  log-normal  model  to  ob- 
tain the  same  probabilities.  The  curves  for  N = 1 and  30  and  Pfa  = 10“6  are  shown  in 
Fig.  20. 

2.  Binary  Integrator  and  Rank  Detector 

The  threshold  values  and  the  P#  can  be  found  with  Eqs.  (44)  and  (45),  by  using 
Eqs.  (63)  and  (64)  for  p{x,  /l).  The  threshold  values  are  given  in  Ref.  45,  and  the  detec- 
tion curves  for  N - 3 and  30  for  the  median  (fe  = jV/2)  and  33d  percentile  (fe  = N/3)  are 
shown  in  Fig.  36.  The  following  observations  can  be  made. 

1.  Whether  the  33d-percentile  value  or  the  50th-percentile  value  (median)  is  better 
depends  on  Pfa  and  Pp.  For  the  values  of  Pfa  and  Pjj  investigated,  there  is  only 
a 1-dB  difference  between  the  two.  Furthermore,  since  the  difference  is  so 
small,  only  a minor  improvement  can  be  made  by  finding  the  optimal  rank. 

2.  For  a fixed  Pq,  as  Pfa  becomes  smaller,  the  33d  percentile  becomes  the  better 
detector. 

3.  For  a fixed  Pfa,  as  PD  becomes  larger,  the  50th  percentile  becomes  the  better 
detector. 

Comparing  the  mean  and  median  detectors  in  Fig.  37,  one  concludes  that  the  mean  re- 
quires a smaller  3/N  than  the  median.  However,  the  difference  decreases  with  increasing 
N.  For  N = 3,  the  mean  is  about  1 dB  better  than  the  median;  for  N = 30,  the  mean  is 
less  than  0.2  dB  better  than  the  median.  The  performance  difference  between  the  two 
detectors  for  large  sample  sizes  can  be  obtained  from  the  ARE.  In  Appendix  C,  it  is 
shown  that  the  ARE  of  the  mean  with  respect  to  the  median  is  1.76.  If  a noncoherent 
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Fig.  37— Probability  of  detection  vs  S/N  for  the  contaminated-normal 
distribution  with  7 = 0.25  and  K = 2.25  [44] 


integration  gain  of  10  log  >/N  is  appropriate  for  small  signals,  the  median  is  approximately 
1.2  dB  better  than  the  mean  for  very  small  signals. 


3.  Trimmed-Mean  Detector 


With  the  previously  outlined  procedure  of  Ref.  49,  a Monte  Carlo  simulation  was 
performed.  For  Ny  = 2 and  N2  = 5,  the  equation 


—log  Pfa  = -1.691  + 0.3123T  + 0.0308T2 , 


(65) 


relating  the  threshold  T to  Pfa,  was  obtained  for  the  contaminated-normal  model,  with 
7 = 0.25  and  K = 2.25.  The  Pjj  curves  were  generated  using  Monte  Carlo  techniques;  the 
results  are  given  in  Fig.  38.  Again,  the  detector  with  the  smaller  ranks  (Ny  - 2 and 
N2  = 5)  is  the  better  detector.  Comparing  the  trimmed-mean  detector’s  performance  with 
that  of  the  mean  (Fig.  35)  indicates  that  the  trimmed-mean  is  the  better  detector,  re- 
quiring 0.5  to  1 dB  less  S/N  than  the  mean. 


4.  Fluctuating  Targets 


The  PD  curves  for  fluctuating  targets  are  generated  by  the  Monte  Carlo  method  dis- 
cussed in  the  previous  section  on  the  log-normal  method.  For  the  contaminated-normal 
model  with  parameters  K and  7,  the  in-phase  and  quadrature-phase  components  are 


58 


VV'.v.v 


> > 


C-  * 


l«4  »•  * K * > 

L ~ - • _ 


% * 
S .S'.  • 


.•.«’.VV.YA 


k 


L.s-i.-vi-s'y 


•V-V-V-*:  •••. 

.;s';-.s. 

r 


i*  -m  ' * * . •-* 

e***  • -5\ 


P>i.v 

• 

-'O 

/■V’* 

/-  i..“ 

is'IvV-* 

»*  . ='■  % 1 

-sswf 

'*2'*->K*J* 

' « -•  %V*  “ »L  y 

> s 


„V 

V vv 


c> 


' ,*  %«**  / V 


- * * * ’ » <*, **  * > 
v v\*V  ■ 


GERARD  V.  TRUNK 


Vi  = [~2in  (w/i )J 1/2  sin  {2mi2)S(ui3)  + A 


z i = [-2  Kti  (ufl)] 1/2  cos  (2irui2)S(uiA) 


S(Uy)  = 


Ka  , uy  < 7 


a , utj  > 7 , 


where  {«,y}  are  independent,  uniformly  distributed,  random  numbers  on  (0, 1).  Using 
Eqs.  (55),  (66)  in  place  of  (06),  (60),  and  (62),  a Monte  Carlo  simulation  was 
performed;  the  results  are  given  in  Ref.  55.  In  Fig.  39,  the  fluctuating  results  are  com- 
pared with  the  nonfluctuating  results.  This  figure  illustrates  the  fact  that  for  a given 
number  of  samples  N,  the  smaller  the  number  of  ranked  samples  involved  in  the  sum 
(n2  - n1  + 1),  the  larger  the  required  increase  in  S/N  to  maintain  the  desired  Pjj.  That 
is,  the  median  suffers  the  largest  “fluctuation  loss,”  the  trimmed  mean  the  next  largest, 
and  the  mean  the  smallest.  Thus,  for  fluctuating  targets  and  contaminated-normal 
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Fig.  39— Comparison  of  various  detectors  for 
Pta  = 10— 6,  for  the  contaminate  Vnormal 
distribution  (7  = 0.25,  K - 2.25),  tnd  for 
fluctuating  and  nonfluctuating  targets  [37] 
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densities,  the  mean  is  the  most  effective  detector,  the  trimmed  mean  is  the  next  best,  and 
the  median  the  least  effective. 


C.  Pulse-to-Pulse  Processing 

Trunk  [39]  investigated  the  problems  associated  with  pulse-to-pulse  processing  of 
data  taken  with  a high-resolution  radar.  His  results  are  summarized  in  this  section.  It 
should  be  noted  that  if  integration  takes  place  during  a single  scan,  the  detection  curves 
(given  in  Secs.  V.A  and  V.B)  cannot  be  used,  because  the  integrated  samples  are  not 
independent  samples  from  p(x).  If  they  are  used,  a higher  number  of  false  alarms  than 
predicted  will  occur.  To  illustrate  this,  the  threshold  values  for  N = 3 and  Pfa  = 10“ 2 , 
and  10"3  were  taken  from  Ref.  10.  Letting  Sj  be  the  sum  of  three  samples  taken  12.5  ms 
apart,  we  calculated  3,000  values  for  Sj  from  the  FHR  data  (run  1153,  HUS  [10]).  The 
sums  were  compared  to  the  thresholds  and  values  of  0.15,  0.016,  and  0.0032  were 
obtained.  Thus,  the  Pfa  is  higher  than  predicted;  the  relative  error  becomes  larger  for 
smaller  Pfa-  The  only  way  of  calculating  the  correct  thresholds  for  a detector  summing 
N pulses  per  scan  is  to  calculate  (approximate,  as  in  Sec.  III.A)  the  density  of  Sj  rather 
than  that  of  x While  Trunk  approximated  a density  of  for  the  HUS  case,  he  noted 
that  adaptive  thresholding  was  a more  fruitful  method. 

A common  adaptive  thresholding  technique  is  the  cell  averaging  constant-ra’ise-a] arm- 
rate  (CFAR)  method  considered  by  Finn  and  Johnson  [57].  Their  detector  for  the  ;'th 
range  cell  can  be  written  as  a ratio, 


i N 

N Z! 


M N 


M L £ (-T?/+8  + l +x?,j-i-l) 


£=1  /-I 

where  Xjj  is  the  ith  envelope-detected  sample  in  the  ;th  range  cell.  (The  samples  on  either 
side  of  tne  test  cell  are  not  used  because  when  a target  is  present,  the  Gaussian  pulse  shape 
causes  target  returns  in  the  adjacent  cells.)  A decision  is  made  by  comparing  Rj  to  a thresh- 
old T(CFAR).  This  detector  provides  CFAR  when  the  noise  samples  are  Rayleigh  distri- 
buted. The  threshold  T(CFAR)  is  a function  of  the  number  of  reference  cells  and  ap- 
proaches the  threshold  T (which  is  the  threshold  if  the  a of  the  Rayleigh  density  known 
a priori)  as  the  product  MN  approaches  infinite.  The  difference  in  S/N  caused  by  2"(CFAR) 
being  larger  than  T is  called  the  “CFAR  loss.”  It  can  be  found  in  Mitchell  and  Walker  [58] 
and  it  is  given  in  Table  14  for  N = 3.  Thus,  to  minimize  the  “CFAR  loss,”  one  wants  M as 
large  as  possible. 

Unfortunately,  if  one  views  the  clutter  in  Fig.  40,  which  was  generated  by  assuming 
the  sea  to  be  frozen  and  converting  time  data  in  Fig.  8 to  spatial  data,  one  sees  that  the 
return  goes  from  a minimum  to  a maximum  in  30  or  40  ft  (9.1-12.2  m).  Therefore, 

Trunk  [39]  suggested  setting  M small  to  detect  in  the  nulls  of  the  sea  clutter  (for  targets 
that  are  not  shadowed,  i.e.,  low-flying  targets)  and  evaluating  the  detector  with  M = 1 and 
actual  clutter  data.  First,  if  x{j  are  independent  samples  from  the  conditional  density 
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P(jc|ct),  which  is  Rayleigh,  then  the  density  of  Rp  under  the  assumption  of  no  signal  present, 
is  the  F-distribution  with  2 N and  AN  degrees  of  freedom.  The  thresholds  can  be  found  in 
Ref.  24.  Since  it  has  been  assumed  that  p(fc|o)  is  a Rayleigh  density,  xfj  has  a chi-squared 
density  with  two  degrees  of  freedom.  If  the  chi-densities  obtained  in  Table  4 were  used, 
the  degrees  of  freedom  for  the  F-distribution  would  be  greater  than  2 N and  AN.  Conse- 
quently, the  threshold  obtained  from  Ref.  24  would  be  lower.  Thus,  we  see  that  the 
Rayleigh  assumption  is  a conservative  approach:  the  actual  Pfa  will  be  lower  than  desired. 


To  calculate  the  Pq  of  the  detector,  a Monte  Carlo  simulation  that  used  actual  clutter 
data  (run  1153,  HUS  [10])  was  performed.  The  output  of  the  cell  containing  signal  was 
generated  by 


Zfj  = xfj  cos0,-  + A)2  + ( xfj  sinfl,^ 


(68) 


where  is  the  actual  FHR  data  sample,  A is  the  nonfluctuating  target  amplitude,  and 
0,-  is  a random  number  uniformly  distributed  on  (0,  2tt).  The  detection  curve  for  N = 3 
and  Pfa  = 10"6  (which  was  derived  from  512  cases)  is  shown  in  Fig.  41.  The  ratio  detec- 
tor is  compared  to  the  one-pulse  detector  (i.e.,  *,•  is  compared  to  an  appropriate  thresh- 
old) and  the  three-pulse  integrator 
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is  compared  to  a threshold,  assuming  that  p(Sy)  is  known  a priori.  (This  is  equivalent  to 
a ratio  detector  with  M = °°).  The  ratio  detector  is  better  than  the  three-pulse  integrator 
except  for  PD  > 0.96,  in  which  case  the  integrator  can  be  1.2  dB  better. 


The  explanation  for  this  is  that  when  signal  strength  is  fairly  low,  the  ratio  detector 
can  detect  signals  in  the  null  of  the  clutter.  On  the  other  hand,  when  signal  strength  is 
high,  resulting  in  high  PD,  signals  must  be  detected  in  all  regions,  including  the  high- 
clutter  region.  However,  in  the  high-clutter  regions  the  three-pulse  integrator  is  better, 
since  no  CFAR  loss  is  suffered  in  estimating  the  threshold.  While  the  PD  is  fairly  low  for 
low  S/N,  it  should  be  noted  that  this  PD  is  for  a single  scan.  It  can  be  improved  by  using 
scan-to-scan  integration  or  a binary  integrator.  It  should  be  noted  from  Table  14  that  the 
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Fig.  41— Comparison  of  the  ratio  detector  with  the  one- 
and  three-pulse  integrators  using  sea  clutter  data  [39] 


63 


I la 


k'  O, 


L 


•;*  \-  ’ • - > -v‘ 

Y VV’-^YY 


v*. 


l 


IVY.-YY -Yaffil 


£■ 


!i  :U- Y*  i;  > 


S 'AvM !}*- 


j Y-  " •*  ”.  *‘m  -*-£4*^ 

Lvi-  v»v‘ 


cv 


wv’-'v-rw 


I 


• •*  V-r.Y5.k- 


% .%  ^ 


GERARD  V.  TRUNK 

ratio  detector  with  M = 1 incurs  a CFAR  loss  of  7.8  dB.  If  the  geometry  of  the  situation 
(the  range  resolution  with  respect  to  the  water  wavelength)  permits  M to  equal  2,  the 
CFAR  loss  can  be  reduced  to  about  4.0  dB. 


VI.  SUMMARY 

When  the  range  resolution  of  a radar  becomes  less  than  the  water  wavelength,  the 
probability  density  of  the  radar  echoes  from  the  sea  has  a longer  tail  than  the  Rayleigh 
density.  This  non- Rayleigh  density  p(x)  has  been  approximated  by  the  log-normal  and 
contaminated-normal  densities.  However,  p(x)  is  neither  of  these  densities.  Rather,  the 
non-Rayleigh  nature  of  p(x)  is  due  to  a spatially  varying  density  p(*|a).  The  density  of 
p(x)  can  be  expressed  as 
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p(x)  = J p(xlo0)p(a0)  do 


where  p(o0)  is  the  probability  density  r.f  a0  and  is  related  to  the  low-frequency  sea  spec- 
trum (i.e.,  large  wave  structure).  The  conditional  density  p(x|a0)  is  a Ricean  density  and 
the  dominant  scatterers  are  associated  with  return  from  breaking,  or  near-breaking,  waves. 
From  the  analysis  of  variance,  the  following  conclusions  were  drawn  about  the  effect  of 
various  parameters  on  the  density  p(x): 

1.  Data  taken  with  horizontal  polarization  have  a large  clutter  spread  than  those 
taken  with  vertical  polarization. 

2.  L-band  data  have  a larger  clutter  spread  than  X-band  data.  This  is  true  for  very 
high  sea  states  and  may  be  true  for  lesser  sea  states. 

3.  Upwind  and  downwind  data  have  a larger  clutter  spread  than  crosswind  data. 

4.  Small-pulse  data  have  a larger  clutter  spread  than  large-pulse  data.  However, 
if  the  pulsewidth  is  smaller  than  the  water  wavelength,  changes  are  no  longer 
significant.  (This  ignores  the  extremely  small  (5  ft  by  5 ft;  1.5  m by  1.5  m) 
resolution  cells.) 

Analysis  of  variance  was  also  applied  to  the  conditional  density  p(x|o0),  and  the  follow- 
ing conclusions  were  reached: 

1.  Data  taken  with  vertical  polarization  follow  more  closely  the  Rayleigh  distribu- 
tion than  horizontal-polarization  data,  which  are  more  Ricean. 

2.  Large-pulse  data  are  more  Rayleigh,  while  short-pulse  data  are  more  Ricean. 

3.  Wind  direction  has  no  effect  on  p(x|a0). 

A simulation  of  the  sea  surface  has  been  developed.  It  has  been  used  to  predict  the 
variation  of  p(x)  and  p(o0)  with  various  radar  parameters.  The  simulation  has  also  been 
used  to  cr’ouiate  the  probability  that  a small  surface  target  is  being  shadowed  and  to 
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show  that  the  times  that  the  target  in  not  shadowed  are  correlated  with  large  clutter 
returns. 
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Various  problems  associated  with  detection  of  targets  have  been  investigated.  Using 
samples  separated  by  about  10  ms  (the  decorrelation  time  usually  stated  for  X-band  sea 
clutter  return)  or  employing  frequency  diversity  produced  independent  samples  from 
p(*|a0),  not  from  p(x).  Independent  samples  from  p(x)  can  be  obtained  only  by  col- 
lecting samples  from  different  scans.  Since  p(x)  has  been  modeled  by  log-normal  and 
contaminated-normal  densities,  det  >ction  curves  for  these  two  densities  have  been  gen- 
erated for  the  mean,  median,  binary  integrator,  and  trimmed-mean  detectors  for  fluc- 
tuating and  nonfluctuating  targets  using  independent  samples  (scan-to-scan  processing). 
Overall,  the  trimmed  mean  detector  is  the  best.  However,  because  of  implementation 
problems,  the  appropriate  binary  integrator  should  be  used.  If  pulse-to-pulse  processing 
is  used,  it  is  recommended  that  if  the  target  is  above  the  surface  a small  number  of  refer- 
ence cells  be  used,  so  that  the  target  can  be  detected  in  the  clutter  nulls.  On  the  other 
hand,  if  the  target  is  a small  target  on  the  surface,  it  can  be  shadowed  and  the  sea  surface 
simulation  must  be  used  to  evaluate  detector  performance. 

The  question  of  what  polarization  should  be  used  remains.  Croney  and  Woroncow 
[59]  recommend  vertical  polarization,  and_  Trunk  and  George  [10]  recommend  horizontal 
polarization.  This  author  believes  the  question  is  now  open;  no  generalization  can  be  made. 
Depending  on  radar  parameters,  data-processing  constraints,  and  environmental  require- 
ments, either  horizontal  or  vertical  polarization  could  be  used.  However,  in  general  this 
author  favors  vertical  polarization  because  of  the  faise-alarm  problem  associated  with  the 
sea  spikes  attributed  mainly  to  horizontal  polarization. 
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Appendix  A 

ANALYSIS  OF  SPATIALLY  VARYING  RAYLEIGH  MODEL 


It  can  be  shown  from  the  frequency-diverse  FHR  data  that  the  conditional  density 
p(x|or,)  is  not  a Rayleigh  density.  If  xt  and  y{-  are  the  independent  backscatter  samples 
from  the  two  frequencies,  the  conditional  densities  of  the  samples  are 


P(*,-|0,)  = ~ exp(-  xf/2af) 
cf 


and 


(Al) 


p(y,K)  " exp(-  yf/2af),  (A2) 

o? 

where  ai  is  a random  variable.  The  same  ai  can  be  used  for  both  frequencies,  since  the 
time  separation  (1/2560  s)  is  very  small  and  the  frequency  difference,  of  the  order  of  the 
reciprocal  of  the  pulsewidth,  changes  oQ  only  slightly. 

Since  a;  changes  rapidly  with  time,  only  a few  samples  of  xi  and  are  available  for 
a fixed  crf.  Consequently,  a ratio  z-  = *./y.,  whose  density  is  easily  shown  to  be 

p(zi)  = 2zi/(z?  + l)2,  (A3) 

is  formed.  The  ratio  p(z{)  is  independent  of  a..  Thus,  even  though  oi  changes,  all  the 
samples  can  be  used  to  test  whether  Eq.  (A3)  gives  the  correct  density  for  the  ratio.  From 
the  FHR  data,  the  sample  distribution  of  independent  (constructed  from  1,024  samples 
taken  12  ms  apart)  was  compared  to  Eq.  (A3)  using  the  Kolmogorov-Smirnov  (KS)  one- 
sample  test.*  The  results  of  the  test  are  shown  in  Table  Al;  D is  the  maximum  difference 
between  the  sample  and  theorized  distributions  and  PD  is  the  probability  that  the  dif- 
ference will  be  less  than  D when  the  theorized  distribution  is  the  true  distribution.  All  12 
cases  are  rejected  at  the  a = 0.002  level  (i.e.,  PD  > 0.998),  indicating  that  xi  and  y;  are 
not  independent  random  variables  of  a conditional  Rayleigh  density.  However,  the  sur- 
prising thing  is  that  the  sample  density  is  more  peaked  (has  a narrower  spread)  than  Eq. 
(A3).  This  indicates  that  either  x;  and  yi  are  correlated,  or  the  density  of  x-  and  yi  is 
more  peaked  than  a Rayleigh  density,  like  a chi  or  Ricean  density. 

However,  it  can  be  shown  that  the  first  of  these  explanations  is  not  possible.  To 
test  whether  correlation  between  je{  and  yt  can  account  for  the  observed  peaked  sample 
density  of  z,  one  must  calculate  the  density  of  z when  x and  y are  correlated  Rayleigh 


*S.S.  Wilks,  Mathematical  Statistics,  Wiley,  New  York,  1962. 
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Table  Al  - Test  of  FHR  Data 
for  a Conditional  Rayleigh 
Distribution 


Identifier 

D 

PD 

VUL 

0.081 

0.999 

HUL 

0.080 

0.999 

VDL 

0.059 

0.998 

HDL 

0.093 

0.999 

HDS 

0.224 

0.999 

VDS 

0.161 

0.999 

VCS 

0.152 

0.999 

HCS 

0.166 

0.999 

VCL 

0.059 

0.998 

HCL 

0.133 

0.999 

HUS 

0.192 

0.999 

vus 

0.136 

0.999  . 

random  variables.  Let  the  in-phase  and  quadrature  components  of  x and  y be  xg,  xc,  y#, 
and  y ; let  the  in-phase  components  be  independent  of  the  quadrature  components  and 
let  both  have  correlation  p.  The  joint  density  of  the  components,  then,  is 


P(xs,  xc , ys,  yc)  = 


exp[-  (x2  - 2px  yc 


.-c  ~r~cJC 

(27r)2(l  - p2) 

+ y2  + X*  - 2 pxv.  + ys)/2(l  - p2)]  . 


(A4) 


Using  transforms  xs  = x sin  0,  xc  = * cos  0,  ys  = y sin  y,  and  yc  = y cos  y,  one  obtains 

1 O O 

- exp{-[x2  +y2 


p(x , y,  0,  y)  = 


(2tt)2(1  - p2) 

- 2pxy  (cos  0 cos  y + sin  0 sin  7)]  /2(1  - p2)} 
Letting  a = 0 - 7 and  (3  = 0 and  integrating  reveal  that 

r2* 

^0 

= f f exP  [2p*y  COS  a/2(l  - p2)]  da  dp 


(A5) 


r27Tr 

J J exp  [2pxy  cos  (0  - 7)/2(l  - p2 )]  da  d7 


= (2t r)2/0 


pxy 


1-P2 
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Table  A2  — Test  of  FHE  Data  for  a 
Correlated  TVR  Distribution 


Identifier 

Optimum  p2 

D 

PD 

VUL 

0.4 

0.015 

0.037 

HUL 

0.3 

0.042 

0.951 

VDL 

0.3 

0.027 

0.558 

HDL 

0.4 

0.037 

0.879 

HDS 

0.4 

0.138 

0.999 

VDS 

0.4 

0.077 

0.999 

VCS 

0.4 

0.069 

0.999 

HC3 

0.4 

0.088 

O.V.*'- 

VCL 

0.3 

0.031 

0.749 

HCL 

0.4 

0.046 

0.974 

HUS 

0.4 

0.104 

0.999 

VUS 

0.4 

0.060 

0.998 

■ 
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0.2  when  the  frequency  difference  is  the  reciprocal  of  the  pulsewidth,*  and  (d)  the 
smallest  “Opt  £2”  was  0.3,  the  hypothesis  that  a:  and  y are  correlated,  Rayleigh-distributed, 
random  variables  must  be  rejected.  It  is  worth  noting  that  the  value  of  p2  was  set  by 
minimizing  D\  a procedure  which  biased  the  test  in  favor  of  acceptance. 
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*V.W.  Pidgeon,  “Time  Frequency  and  Spatial  Correlation  of  Radar  Sea  Return,”  Proc.  Amer.  Astronaut. 
Soc.  Symp.  (Boston,  May  25-27,  1967),  in  Use  of  Space  Systems  for  Planetary  Geology  and  Geophysics 
The  Society,  Tarzana,  Calif.,  1968,  pp.  455*458. 
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Appendix  B 

BRIEF  DESCRIPTION  OF  ANALYSIS  OF  VARIANCE 


The  analysis  of  variance  (ANOVA)  procedure  will  be  demonstrated  by  finding  the 
variation  of  the  clutter  distribution  with  changing  frequency  and  azimuthal  squint  angle. 
As  shown  in  Table  Bl,  the  distribution  can  be  represented  by  a single  number  o.  In 
analysis  of  variance,  the  data  o{jk  are  represented  by  a linear  model  which  consists  of  a 
mean  n,  a frequency  effect  F,  an  azimuthal  angle  effect  A-,  and  a random  error  e,./; , such 
that 


i = 1, ...,  Np 

aijk  = v + Fi+Aj+  eijk  ; = 

k=l,...,NR 


where  Np  is  the  number  of  frequencies,  NA  is  the  number  of  azimuthal  angles,  NR  is  the 
number  of  repetitions  of  each  frequency  azimuth  ^ jse, 


Nr 


f F,  = 0,  £ A,- 0. 

i=l  7=1 


Na 


and  el]k  are  independent  Gaussian  random  variables  with  mean  0 and  unknown  variance  X2. 
The  significance  of  a change  in  frequency  is  checked  by  tes  ing  the  hypotheses, 


Hq  :Fi  = 0,  for  all  i (no  frequency  effect  present) 
H1  :F;  =£  0,  for  some  i (frequency  effect  present), 


using  an  F-test  (an  optimal  test  for  the  equality  of  unknown  variances).  The  test  is  per- 
formed by  taking  the  ratio  R of  two  statistics;  Sp,  an  unbiased  estimator  of  X2  + 4(FX  - 
F2)2,  and  S|,  an  unbiased  estimator  of  X2.  Under  HQ,  the  ratio  has  an  F-distribution, 

J 1 TT  X*  1 * x 1 n i 1 i; n _ i r - A •*  n _ 


and  under  Hy  the  ratio  has  a noncentral  F-distribution.  For  a type-/  error  of  0.1,  the 


threshold  value  is  found  in  Duncan  to  be  3.2.*  In  Table  B2,  R = 55.8;  consequently,  the 
null  hypothesis. is  rejected.  That  is,  a frequency  effect  is  present.  The  procedure  is  re- 
peated for  the  azimuthal  angle,  and  the  effect  is  not  significant.  Interactionsf  were  found 
not  to  be  significant  at  the  0.1  level  in  this  or  any  other  data  used  in  this  paper.  Conse- 
quently, to  avoid  unnecessary  complications,  the  results  of  interactions  were  not  men- 
tioned in  the  report.  Further  details  about  analysis  of  variance  can  be  found  in  either 
Duncan*  or  Fisher.-*-. 


♦A.J.  Duncan,  Quality  Control  and  Industrial  Statistics,  Irwin,  Homewood,  DL,  1959. 

'R.A.  Fisher,  Statistical  Methods  for  Research  Workers  Oliver  and  Boyd,  Edinburgh,  1941. 
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4 FR  Radar 
Parameters 


Azimuthal  £ 

quint  Angle 

0° 

15° 

CO 

o 

o 

45° 

Frequency 
L band 
X band 


f Analysis  of  Variance:  The  Effect  of  Various 
rs  on  the  4FR  Clutter  Distributions* 


Effect 

Mean  Square 

Variance  Ratio 

Conclusion 

Frequency 
Azimuthal  angle 
Sampling  error 

S|  = 5.640 
Sf  = 0.062 

sf  = 0.101 

Sj/Sl  = 55.84 
Sl/Sl  = 0.61 

Frequency  effect  is  present 
Azimuthal  effect  not  present 

♦From  reference  [16]. 
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Appendix  C 

ASYMPTOTIC  RELATIVE  EFFICIENCY  OF  MEAN 
AND  MEDIAN  DETECTORS 


The  asymptotic  relative  efficiency  (ARE)  can  usually  be  calculated  rather  easily  by 
employing  a concept  known  as  efficacy.*  That  is,  the  ARE  of  two  detectors  and  d2 


is 


ARE  (d2>  cfj)  = 


§2 


(Cl) 


where  S.  is  the  efficacy  of  detector  d..  Specifically,  given  a binary  hypothesis-testing 
problem  (HniA  = A0  vs  H^.A  > A0),  if  the  detector  df  is  based  on  a statistic  Tin  - T,(xi,x2, 
...,  xn)  and  if  E{T.ni  = \pin (A)  and  Var{2]n)  = a?  (A)  and  m and  8 are  defined  by 


where 


= *f?,M  o)  > 0 


(C2) 


*!?>Wo)s 


3A 


m 


and 


lim  n"m5 

n-*°° 


*6°  Wo> 


= C>0, 


(C3) 


then  the  efficacy  of  detector  d,-  is 


6 = 


Il/m5 


(C4) 


subject  to  certain  regularity  conditions.'*’  When  dealing  with  the  class  of  translation 
alternatives,  the  following  conditions  generally  apply:  m = 1,  8 = 1/2,  and  the  evaluation 
of  Eq.  (C4)  is  very  simple.  However,  when  dealing  with  other  alternative  classes  (e.g., 
noncoherent  detection),  generally  m = 2,  8 = 1/4,  and  the  evaluation  of  Eq.  (C4)  is  not 
always  simple. 


♦E.J.G.  Pitman,  “Lecture  Notes  on  Nonparametric  Statistical  Interference,”  Columbia  University,  Sprinj 
19-48. 

'G.E.  Noether,  “On  a Theorem  of  Pitman,”  Ann.  Math.  Statist  26,  64-68  (1966). 
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An  expression  for  the  efficacy  of  the  mean  against  the  log-normal  density  is  derived 


u jl  r 

9 A2  *4 


xp(x,  A)  dx 


by  taking  the  second  derivative  inside  the  integral  and  dropping  terms  that  when  evaluated 
at  zero  signal  strength  are  zero.  The  efficacy  of  the  mean  can  be  written  as 


j4(o2-l)  . 16(a2  fin*  + (fin  x2) 

W — 7T~  + 71 


X (27ra2)-1^2exp  [-  2 (fin  x)2/o2]  dx  , 


„a2/„a2 


where,  for  the  efficacy  of  the  sample  mean,  F,  (x)  = nx,  T1  = °°,  and  al  - [nea  (e  - 

1)]1/2. 

The  expression  for  the  efficacy  of  the  median  is  derived  by  first  noting  that  making 
a decision  when  the  median  is  greater  than  a threshold  value  T is  entirely  equivalent  to 
making  the  decision  by  counting  the  number  of  samples  greater  than  T.  Thus,  as  the 
number  of  samples  approaches  infinity,  T approaches  the  median  value  of  the  distribution, 
since  the  median  detector  is  consistent.  Hence,  if 


rxm 

E{Tn}  -np  = n J p(x,  A)  dx, 


the  efficacy  of  the  median  is 


_ „ r<aw>f 

-n  , 

. 0.5 y/n  . 


which  reduces  to  Eq.  (C6)  with  F2(R)  = 1,  T„  = 1,  and  o2  = 0.5A/7T 


Now,  if  one  repeats  the  previous  procedure,  one  obtains  for  the  efficacy  of  the  mean 
against  the  contaminated-normal  density 


&r[JQ 


tt-7)2 


_ 1 + exp 


+ -2^-  (- 1 + — ] exp  (——— 
K4o*  \ 2 K2o2  \2 K2o2 


(C9) 
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27(1  - 7) 
Kiro 4 


n- 


i + 


x2  cos2  0 


X exp 


+ JT1  -1  + 


X exp 


- [x2  + (K2-l)x2  cos2  0] 

2 K2o2 
x2  cos2  0 
K2o2 

~(x2K2  - {K2  - l)x2  cos2  0j\ 
2K2<j2 


where  Fjfc)  = nse,  Tj  - and  Oj  is  the  standard  deviation  of  the  sample  mean  for  the 
noncoherent  contaminated-normal  distribution  with  A = 0,  which  is  (12  - ir/2)1^2  <j\/n 
for  7-O  and  is  1.11  a\/n  for  7 = 0.25  and  K = 2.25.  The  latter  value  for  was 
found  by  numerical  integration. 

The  efficacy  of  the  median  against  the  contaminated-normal  density  reduces  to  Eq. 
(C9)  with  *2{x)  = nx,  T2  = the  median  value  of  the  noncoherent  contaminated  normal 
distribution  (which  Is  (2  in  2)^2  0 for  7 = 0 and  is  1.41o  for  7 = 0.25  and  K = 2.25), 
and  a2  = 0.5 y/tt. 


For  the  Rayleigh  case  (7  = 0),  the  integral  in  Eq.  (C9)  can  be  performed.  It  yields 

= 0.91 5n 


1 

8 m 

&1  = 2(1  - it  14) 


and 


6-  = (fin  2)2n  = 0.48n 


for  the  efficacies  of  the  mean  and  median,  respectively.  The  ARE  is 

&2 

ARE  (median,  mean)  = — = 0.525. 

sl 

This  implies  that  in  the  limit  as  A -*■  0,  in  order  to  maintain  the  same  probabilities  of 
false  alarm  and  detection,  the  mean  requires  only  52.5  percent  as  many  samples  as  the 
median.  Thus,  the  mean  is  the  better  large-sample  detector  for  the  Rayleigh  distribution. 

For  the  noncoherent  contaminated-normal  (7  0)  and  the  log-normal  densities, 

Eqs.  (C6)  and  (C9)  have  been  evaluated  numerically;  the  results  obtained  are  given  in 
Table  Cl.  These  show  that,  in  all  cases,  the  median  is  the  better  large-sample  detector. 
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Table  Cl  — Asymptotic  Relative  Efficency  of  the  Mean 
and  Median  Detectors 


Distribution 

Efficacy  of 
Sampled  Mean 

Efficacy  of 
Median 

ARE 

Log-normal  (a  = 3) 

0.286n 

5.334n 

18.6 

Log-normal  ( a = 6) 

0.0103n 

1.333n 

129 

Contaminated-normal 

(7  = 0.25  and  K=  2.25) 

0.0137n 

0.0238n 

1.76 

E.  K.  AL  Hussein*  has  solved  the  integral  in  Eq.  (C9)  for  S2  and  has  )btained 

S2  = S/i/tto2.  (CIO) 

Substituting  o = 0.2303o(dB)  into  Eq.  (CIO)  with  a = 3 dB  and  a = 6 dB,  one  obtains 
5.334n  and  1.333/1,  respectively. 


Appendix  D 

IMPORTANCE  SAMPLING 


A straightforward  method  of  determining  the  required  threshold  for  a given  Pfa  is  to 
perform  a Monte  Carlo  simulation.  Unfortunately,  for  P.  - 10"6,  more  than  one  million 
repetitions  would  need  to  be  run,  and  the  computation  time  on  any  computer  would  be 
large.  However,  a simulation  that  uses  importance  sampling  can  be  used.*  The  main  pur- 
pose of  importance  sampling  is  to  modify  the  probabilities  that  govern  the  outcome  of 
the  basic  experiment,  of  the  simulation  so  that  the  event  of  interest  (i.e.,  a false  alarm) 
occurs  more  frequently.  This  distortion  is  then  compensated  for  by  weighting  each  event 
by  the  ratio  of  its  probability  if  the  true  probabilities  had  been  used  to  its  probability 
with  the  distorted  probabilities.  Consequently,  by  proper  choice  of  distorted  probabili- 
ties, the  number  of  repetitions  can  be  reduced  greatly.  For  instance,  in  estimating  the 
mean  of  a function  Q(x) 


£[&(*)]  = fax)  dP(x), 


(Dl) 


where  P(x)  is  the  distribution  of  x;  the  mean  can  be  estimated  by  selecting  m independent 
samples  x-  from  P(x)  and  associating  the  probability  1/M  with  each  event.  Then,  i?[Q(x)] 
can  be  estimated  by 


1 M 

u L «*()■  <D2) 

»=1 

The  importance-sampling  technique  used  the  Radon-Nikodyn  derivative  to  express  the 
mean  value  of  Q(x)  by 


TOW! 


(D3) 


where  G(x)  is  a distribution  function.  Now,  Ji[Q(x)]  can  be  estimated  by  selecting  M 
independent  samples  from  G(x)  and  associating  the  probability  dP(xi)/MaG(xi ) with 
event  Q(x{).  Thus,  F[Q(x)]  is  estimated  by 


i ™ dP{x{) 

if  A QiM  ,?G(x.)  ' 


1=1 


(D4) 


*F.S.  Hillier  and  G.J.  Lieberman,  Introduction  to  Operation  Research  Holden-Day,  San  Francisco,  1967, 
pp.  457-469. 
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Since  Eqs.  (D3)  and  (D4)  are  both  unbiased  estimates  of  Q(x),  it  is  possible  to  select 
G(x)  so  that  the  variance  of  Eq.  (D4)  is  less  than  the  variance  of  Eq.  (D2). 


To  illustrate,  let  us  generate  the  distribution  function  (and  hence  curve  of  threshold 

T ’ 

If 

such  that 


xi  - exp{a[-  2 Cn^-j)]1^2  sin  (27ra|2)}» 


(D5) 


have  a log-normal  density,  with  xm  = 1.  Recalling  that  the  0 m Eq.  (D5)  is  m natural 
units  (0  = 0.2303a  (dB)),  choose  oH,  which  is  greater  Man  0.  To  simulate  die  trimmed 
mean  for  n.  = 2,  n2  - 5,  and  N = 10,  generate  10  samples  xp  usmg  Eq.  (D5)  with  a 
replaced  by  aw  Order  the  fr.}  and  form  the  sum 


"2 

E *<• 

i-n. 


Since  this  procedure  is  repeated  M times,  denote  log-normal  samples  by  x--  and  trimmed- 
mean  samples  by  S-.  The  estimated  cumulative  distribution  of  S-  for  the  log-normal 
model  with  parameter  0 is 


M 


1 m 

5 E Vv 

y=i 


(D6) 


where 


8r 


(1  Sj<T 


(D7) 


(0  Sj>T 


N 2 exp  {-  2[8n(xf.,)]2/o2}/(27ra2xf;.)1/2 

Pj-U  — 


Jm\  2 exp {-  2[in(xij))2lafI}l(2naflxij)112 


This  reduces  to 


(D8) 


While,  a priori,  the  desired  value  of  oH  ^unknown,  an  appropriate  value  can  easily  be 
found  since  the  variance  of  the  estimate  p(T)  is  given  by 


*S.F.  George,  “The  Detection  of  Nonfluctuating  Targets  in  Log-normal  Clutter,”  NRL  Report  6796,  Oct. 
4, 1968. 
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